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ABSTRACT 


The  evaluation  of  the  direction-finding  capabilities  of  non-uniform 
arrays  is  approached  via  algorithm-independent  lower  bounds  on  achievable 
angle  estiaatlon  errors.  Two  classes  of  bounds  are  considered.  The  first* 
known  as  the  Craaer-Rao  bound,  applies  only  to  unbiased  estimates.  Coepact 
analytical  egressions  for  these  bounds  are  developed  which  are  applicable  to 
very  general  direction-finding  probleas,  including  an  arbitrary  nuaber  of 
ealtters. 

It  is  well-known  that  Craaer-Rao  bounds  are  overly  optlnistic  at  low 
signal-to-nolse  ratios.  As  this  ratio  decreases,  a  point  Is  reached  at  which 
estiaatlon  accuracy  decreases  abruptly.  Another  class  of  bounds,  known  as 
Zlv-Zakal  bounds,  provide  lnforaatlon  about  the  location  of  this  threshold 
point.  Their  study  suggests  that  poor  direction-finding  perforaance  occurs  in 
situations  where  the  ealtter  direction  vectors  are  part  of  a  set  which  is 
nearly  linearly  dependent. 

Such  linear  dependence  does  not  occur  in  the  case  of  uniformly  spaced 
linear  arrays  (without  getting  lobes).  However,  It  does  occur  when  eleaents 


are  reaoved  froa  such  arrays.  A  systeaatlc  test  Is  developed  to  test  a  given 
array  geometry  for  this  condition.  -7 

Finally,  direction-finding  perforaapce  of  the  QUICK  LOOK  array  is  evalu¬ 
ated  via  both  bounds  and  simulation.  / 
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I.  INTRODUCTION 


The  evaluation  of  the  direction-finding  capabilities  of  non-uniform  ar¬ 
rays  Is  approached  via  algorithe-independent  lower  bounds  on  achievable  angle 
estimation  errors.  TWo  classes  of  bounds  are  considered.  The  first,  known  as 
the  Craaar-Bao  bound,  applies  only  to  unbiased  estimates.  Since  almost  all 
"reasonable"  estimates  are  asymptotically  unbiased  as  the  slgnal-to-noise  ra¬ 
tio  becoaes  large,  these  bounds  are  useful  for  evaluating  performance  in  that 
situation.  Coapact  analytical  egressions  for  these  bounds  are  developed 
which  are  applicable  to  very  general  direction-finding  probleaa.  Including  an 
arbitrary  number  of  ealtters. 

It  Is  well-known  that  Craaer-Kao  bounds  are  overly  optimistic  at  low 
signal-to-noise  ratios.  As  this  ratio  decreases,  a  point  is  reached  at  which 
estimation  accuracy  decreases  abruptly.  Another  class  of  bounds,  known  as 
Zlv-Zakal  bounds,  provide  information  about  the  location  of  this  threshold 
point.  Calculation  of  these  bounds  requires  numerical  maximisation  of  a  func¬ 
tion.  In  the  case  of  more  than  one  emitter,  this  function  1b  multi¬ 
dimensional  and  the  computational  load  appears  to  be  prohibitive.  Despite 
this  difficulty,  the  study  of  Zlv-Zakal  bounds  suggests  that  poor  threshold 
behavior  occurs  In  situations  where  the  emitter  direction  vectors  are  part  of 
a  sat  which  Is  nearly  linearly  dependent. 

The  simplest  example  of  linear  dependence  Is  a  grating  lobe  of  a  uni¬ 
formly  spaced  linear  array.  When  the  array  elements  are  separated  by  more 
than  half  a  wavelength  (1/2),  there  exist  two  (or  more)  directions  which  pro¬ 
duce  the  same  direction  vector.  It  is  Impossible  to  determine  from  the  array 
data  whether  the  signal  is  from  direction  A,  direction  8,  or  both.  When  the 
spacing  Is  less  than  1/2,  it  can  be  shown  that  for  every  set  of  distinct  di¬ 
rections,  the  corresponding  set  of  direction  vectors  is  linearly  independent, 
and  therefore  no  ambiguity  problems  exist. 

A  special  case  of  non-unlformly  spaced  linear  arrays,  called  thinned  lin¬ 
ear  arrays,  are  those  obtelned  by  removing  elements  from  a  uniform  linear  ar¬ 
ray.  These  exhibit  ambiguities  of  a  more  complicated  nature.  Tor  example, 
when  a  sat  of  three  direction  vectors  Is  linearly  dependent,  the  observed  data 


could  have  cow  froa  direction  A  end  B,  or  A  and  C,  or  B  and  C ,  or  all  three. 
A  given  array  geonetry  can  be  systematically  tested  £or  the  presence  of  line¬ 
arly  dependent  direction  vectors. 

Finally,  direction-finding  performance  of  the  QUICK  LOOK  array  is  evalu¬ 
ated  via  both  bounds  and  simulation.  Results  show  that  the  currently  used 
"phase  only"  processing  is  essentially  optimum  for  a  single  emitter.  For  two 
emitters,  a  processing  algorithm  called  MUSIC  provides  nearly  optimum  perform¬ 
ance  except  in  certain  difficult  cases.  Further  work  is  needed  to  determine 
whether  these  failures  are  due  to  the  array  geometry  or  to  a  weakness  in  the 
algorithm  Itself. 


II.  PERFORMANCE  BOUNDS 

The  first  question  that  erleee  la  the  study  of  unequally  spaced  linear 
arrays  Is  that  of  array  geometry.  What  sort  of  array  geonatrles  are  "beat" 
for  locating  nultlple  ealtters?  Evaluation  of  an  array  via  simulation 
requires  the  choice  of  an  algorithm;  If  performance  Is  bad.  It  nay  be  due  to  a 
poor  choice  of  algorithm.  The  study  of  algorithm  -  Independent  performance 
bounds  Is  motivated  by  the  desire  to  separate  these  two  factors. 

The  oldest  and  best  known  bound  Is  the  Craner-Rao  bound.  It  Is  valid 
only  for  unbiased  estlaators  (or  estimators  of  known  bias).  When  the 
observations  from  which  an  estimate  Is  to  be  made  are  of  the  form 

£. "  A»  <«  >  ♦  1 

where  ^  is  a  deterministic  vector  which  depends  on  a  set  of  unknown  parameters 
a,,  and  ^  is  a  Gaussian  noise  vector,  It  can  be  shown  that  there  exists  an 
estimator  which  achieves  the  Crane r-Rao  bound  asymptotically  as  the  slgnal-to- 
nolse  ratio  A2/e|ca|2  becomes  large. 

As  slgnal-to-noise  ratio  decreases,  a  point  Is  reached  at  idilch  the 
accuracy  of  the  estimate  degrades  abruptly.  This  value  of  S/N  Is  referred  to 
as  the  threshold  S/N;  In  general.  It  will  depend  on  the  estimation 
algorithm.  The  Crame r-Rao  bound  gives  no  Information  about  this  behavior.  A 
bound  which  does  give  such  Information  is  the  Zlv-Zakal  bound. 

1.  The  Craner-Rao  Bound 
A.  General  Discussion 

Derivations  of  the  multlparaneter  Craner-Rao  bound  are  readily  available 
in  the  literature  (e.g.,  [1]).  The  result  may  be  stated  as  follows:  let  a,  be 
a  real  K-dlnensional  vector  parameter,  and  let  p(,s/a)  be  a  probabilistic 
mapping  from  this  parametsr  space  Into  an  observation  space  whose  elements  are 
N-dlmanslonal  real  or  complex  vectors  ».  The  Crams r-Rao  bound  states  that  the 
covariance  matrix  of  any  estimate  a  of  a  satisfies 


9a  T  .  9a 
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i.e.  the  difference  of  theee  two  aetrlcee  Is  positive  seal definite.  Here 


ii  ■  E(a) 
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is  called  the  Fisher  information  aatriz. 

9a 

For  unbiased  estlaates,  —  ■  I  ,  the  identity 

“  oCt_ 

becoaes  siaply 


aatriz,  and  the  bound 


Cov  (i)  >  F_1 


(2.2) 


The  remainder  of  this  section  will  discuss  only  unbiased  estimates. 


It  Is  frequently  the  case  that  only  soae  of  the  unknown  parameters, 
denoted  by  6,  are  of  interest.  The  remaining  paraaeters  y ,  celled  "nuisance" 
parameters,  are  not  of  Interest  but  aust  be  estlaated.  The  Information  aatriz 
can  then  be  partitioned  as  follows: 
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It  Is  readily  shown  that 
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(2.3) 


Ezpresslons  for  the  other  subastrlces  are  easily  obtained,  but  will  not  be 

*  —1 

required  here.  Since  Cov(a)>F  ,  it  follows  that 


•"tiWV  Vn  ?T«> 
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(2.4) 
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for  otherwise  a  principal  alnor  of  Cov  (a)  -F_1  would  be  negative. 

The  aatrlx  whose  Inverse  appears  on  the  right  side  of  (2.4)  is  often 
called  the  reduced  information  aatrlx  for  the  parameter  set  6 ,  denoted  by 

It  consists  of  the  appropriate  subaatrix  of  F  reduced  by  a  positive 
semldeflnite  matrix. 

Application  of  this  result  to  the  case  of  a  single  parameter  shows  that 
the  variance  in  the  estimate  of  a  single  parameter  must  exceed  the  Inverse  of 
the  corresponding  diagonal  element  of  F,  i.e., 

Cov(6)  >  -±-  (2.5) 

60 

B.  Estimation  of  Signal  Parameters  in  Gaussian  Noise 

An  important  special  case  to  which  the  Cramer-Rao  bound  may  be  applied  is 
the  estimation  of  the  parameters  of  a  deterministic  signal  observed  in  the 
presence  of  additive  Gaussian  noise,  i.e., 

z_  ■  ±(a)  + 

where  £  is  a  circular  complex  Gaussian  vector  with  zero  mean  and  covariance  A . 
The  logarithm  of  the  likelihood  function  is  then  of  the  form 

In  p(z/a)  — Mfctm-Jln  jA  |-(£  -  s(a))H  A  *(z-g(a)) 

with  derivative 
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The  Fisher  information  matrix  is  then 


F 
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'  s  H 

HI*? 


(2.6) 
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It  Is  often  the  case  that  some  of  the  real  paraaeters  a_  are  actually  the 
real  and  Imaginary  parts  of  comp lex- valued  paraaeters,  e.g.  unknown  coaplex 
amplitude  of  the  signal  vector  b_.  Denote  this  set  of  coaplex  paraaeters  by 
the  vector  c^  and  assume  that  j»(c)  Is  an  analytic  vector  function  of  c_  i.e., 
each  coaponent  of  _s  is  an  analytic  function  of  each  of  the  coaponents  of  c. 


Define  C 


r^j= 

9c 


It  follows  from  the  Cauchy-Riemann  equations  that 


r*_L-c.  r'hJZ 


Let  now  denote  the  remaining  real  paraaeters  and  define  A“  —  .  The  Fisher 

T  T  \  f  A 

information  matrix  for  the  (transposed)  parameter  vector  (j»  ,  Re  £  ,  Im  _c  ' 


F-  2  Re 


AHA  AHC  jAHC 
CHC  jCHC 

n 

Hermit lan  C  C 


The  reduced  information  matrix  for  the  parameter  vector  alone  Is  then. 


Re  CHC  -Im  CHC  /Re  CHa\ 
Im  (^C  Re  CHC  yim  CHA/ 


F(_)  -  2{Re  A«A  -  (Re  AHC,  -  ImA® 


Evaluation  of  this  expression  Is  greatly  simplified  by  making  use  of  the 
Isomorphism  between  multiplication  of  complex  NxM  matrices  and  of  real  2Nx2N 
matrices  of  special  form. 


tRe  A  -Im  a1  fRe  B  -Im  b!  fRe  AB  -Im  AB 

Im  A  Re  A  In  B  Re  B  In  AB  Re  AB 

A  B  J  -  L  AB 


*We  use  the  letter  J  to  denote  /-I  throughout. 


4# 


It  follows  immediately  that 


‘Re  B  -Ie  b|  _1  ["r«(b”1) 

la  B  Re  B  L^b”1) 


-I«(B_1) 

Re(B_1)_ 


and  that 


(Re  A,  -  IaA)  [fie  B  -la  b]”1  [Re  c]  -  Re(AB_1C) 

[la  B  ReB  J  [la  Cj 

Application  of  these  results  to  the  above  expression  yields 

F*“*-  Re  { AHt I-C(  C^CJ “1 C® J  A}  (2.7) 

and  the  Craaer-Rao  bound  on  the  covariance  aatrlx  of  any  unbiased  estlaate 

A 

u  of  u  is 


C°v  («J  >  t O”1 


(2.8) 


C.  Multiple  Bleaeatary  Signal*  in  Gaussian  Noise 

This  Is  actually  a  special  case  of  the  preceding  section.  The 
observation  Is  modelled  as 

K 

*  - );  ♦  £  *  v(«>  £ * 1 
k»l  * 

where 


,TC.)  . J2*“SA, 

ST 

is  a  noraallsed  direction  vector  which  represents  the  response  of  an  Ideal 
linear  array  with  eleaents  at  positions  xa  to  a  plane  wave  from  direction  u. 
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It  1*  convenient  to  define  the  origin  of  the  x  coordinate  eo  that  £  x  •<)  . 
The  "eleaentary"  signals  are  all  of  the  fora  _v(u)  and  have  unknown  coaplex 
amplitudes.  The  results  of  the  last  section  are  directly  applicable  with  the 


substitution 


Thus, 


where 


£  (a)  -  V(u>2_ 


A  “A~ 1/2  A’1/2^  w 

C  -A~  ^2  3~  .  2y 


3V  3-  3  v 

^  3u  ^  3uj  3uj  ***  Su^  ^ 


.  [Pl-  0 

L  **} 


The  reduced  Infornatloa  watrlx  for  u  Is  then  obtained  using  (2.4).  The  result 


where 


ri'5  -  2  Re  |[£jH  (W  -  ft®W_1ft)  [£]} 


4  V®  a"1v 


ft  4  v®  a'H 


(2.9) 


(2.10a) 


(2.10b) 


(2.10c) 


The  &  notation  la  mnemonic  only  and  ahoold  not  ba  confused  with  that  for 
derivatives. 

The  Craaer-Rao  bound  on  tha  covarlaoca  aatrlx  of  any  unblaaad  aatlnata  of 
ti  la  than  given  by 

Cov(u)  >  j  Re"1  U^jYw  -  ft8*-1*)!*)}  (2.11) 


D.  Extanalon  to  Multlpla  Obaarvatloua 

In  tha  nultlpla  look  eaae,  tha  quaatlon  of  algnal  eoharanca  arlaea.  Two 
Uniting  caaaa  ara  of  particular  Intaraat.  In  tha  Incoharunt  case,  the  algnal 
anplltudas  are  conpletaly  Independent  fron  one  look  to  tha  next;  they  are 
treated  as  additional  unknown  parameters.  The  nth  observation  la  modelled  as 

*.  •  «S>£,  +  i. 

The  most  concise  derivation  of  tha  Information  aatrlx  for  this  case  makes  use 
of  the  Kronecker  product  [2]  of  two  matrices,  defined  by 


•llB  *12B 

*1# 

*21B 

■mib 

V 

If  A  la  N  X  N  and  B  la  R  X  S,  their  Kronecker  product  la  MR  X  MS.  The 
following  easily  verified  properties  of  tha  Kronecker  product  will  be 
required. 


(»B)H  -  lW 


(MB)  (CRD)  -  ACRBD 


If  A,C  and  B,D  are  conformable 


CAM)"1  -  A-1*!’1 


Following  the  development  of  the  preceding  eectlon 
e^Co)  -  VCu)^ 


3J« 

A  «  _Z2£ 

\  du 


6  A 


-1/2 
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Define  the  concatenated  vectore 


Now 


4  -  <4-  4-  -4) 
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where 
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end 
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A  "1/2V  u 

A  V*V 

o 

• 

• 

i _ 

V^v 

-  m"l^V  (2.13) 


From  (2.12)  and  (2.13)  there  follows,  again  using  tha  notation  of  (2.10) 


AHA  -  QH(ItA~  1/2  V)H  (I*A~  ^  V)Q  -  QH(I*Sr)Q  (2.14a) 

C^A  -  (It Id”  1/2  V®)  (ItA~  ^  V)Q  -  (Ifl&)Q  (2.14b) 

(t^C)  -  (ItlT^V)®  dtA-1^V)  -  (Itw)  (2.14c) 

Us«  of  (2.14)  in  (2.7)  gives 

2  Ref  QHf  MW  -  (  »&)H  (  »W_1)  ( I»W)  J  q} 

-  2  M  l  W"1  *0  ll>J}  <2-15) 


which  la  the  reduced  Information  net r lx  for  the  multiple-look  Incoherent  case. 

In  the  other  Halting  case  of  coapleta  coherence,  the  amplitude  and  phase 
relationships  among  the  signals  are  fixed  from  look  to  look;  only  the  overall 
complex  amplitude  changes.  Thus  the  n**1  vector  observation  is  modeled  as 


**  "  V(«)£  +  £  n  n-l.Mjhj-l 

The  number  of  unknown  parameters  in  this  model  is  2(R-1)+R+2JC,  as  opposed  to 
2NR4R  in  the  incoherent  case. 


Proceeding  as  before 

•*(«)  "  V(«).£ 


11 


1 


*a*  -  io  •WVeI  -  Ali*  "til 


when 


there  follow 


lwedletely 


AHC(  c®c)_1cha  -  AHC  (c^c)"1 


»l£j 


-  bBj  Jjj"1]  &t £l 

- 


end  finally,  froa  (2.7) 
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f^)  -  2  d2u{l£jH(w  -  (2.16) 

which  Is  Che  reduced  Inf  one  Cl  on  matrix  for  the  Multiple- look  coherent  cese. 
The  Cramer-Rao  bound  on  the  error  In  estimating  u_  le  again  the  Inverse  of  this 
matrix. 

Although  the  derivation  In  the  two  cases  Is  quite  different.  It  Is  a 
curious  fact  that  (2.16)  can  be  obtained  formally  simply  by  waking  the 
substitution  p|1*bnp  In  (2.15). 

The  principle  results  of  this  section  are  equations  (2.15)  and  (2.16), 
which  are  the  reduced  Information  matrices  for  estimating  the  directions  of 
arrival  of  multiple  emitters  In  interference  of  an  arbitrary  nature,  from 
multiple  independent  observations.  Equation  (2.15)  is  for  the  Incoherent  case 
In  which  the  unknown  complex  amplitudes  of  the  emitters  change  Independently 
from  observation  to  observation.  Equation  (2.16)  covers  the  completely 
coherent  case  in  which  the  unknown  complex  amplitudes  have  the  same  fixed 
relationship  on  every  observation;  the  only  fluctuation  Is  a  complex  scale 
factor  which  changes  from  one  observation  to  the  next. 

E.  Examples 

Single  Emitter  in  White  Holes 

For  a  single  emitter  In  white  noise  with  covariance  o2I,  equations  (2.10) 

2vx 

become  simply  W-1.&-0,  W**~  \  ( — ~)2  -  (2w£)2  .  The  Creme r-Rao  bound 

(2.11)  Is  then  simply 

2  1 

Cov  (u)  >  - jjr — x - i— ,  (2.17) 

2  "  W 

The  covariance  of  u  la  seen  to  be  Inversely  proportional  to  slgnal-to-nolse 
ratio.  Independent  of  the  true  direction  u,  end  dependent  on  array  geometry 
only  through  the  "RMS  array  length"  1. 

Two  Emitters  In  White  Holes 

Bounde  for  the  case  of  two  emitters  with  directions  uj,  u2  are  more 
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complicated.  Kaaulta  have  bean  obtained . for  a  aequence  of  caaea  of  increasing 
complexity.  The  amplitude ,  phase,  and  direction  of  the  flrat  emitter  ere 
aaeumed  to  be  unknown  parameters.  The  corresponding  parameters  for  the  second 
interfering  emitter  are  known  or  unknown,  as  specified. 


(1)  Interference  Known 

If  all  of  the  parameters  of  the  second  emitter  are 
known,  it  can  be  subtracted  out  and  does  not  degrade 
the  accuracy.  Thus  we  get  the  single  emitter  bound. 


Cov  (u,  >  >  - - *■  — - - 

2  l  IpJ2  (2wl) 


2  al 


(2)  Interference  Amplitude  end  Phase  Unknown 


Cov(Uj)  >  Oj 


ivSTT 

'—I  -s' 


(2.18) 


(2»i)  (1-  |v“  vj  ) 

The  function  n(us,  Uj)  is  an  efficiency  factor  (OCn<l)  which  depends  on 
the  source  directions  and  the  array  geometry. 

(3)  Interference  Amplitude,  Phase,  and  Direction  Unknown 


Cov(u.)  >  - 


4  2  1 


H*  x2 


— ij- 

(2ri)2n  ~1_I  1  -  |vj  vx|2 


1  n  l-leZ{ye^} 


(2.19) 


The  last  factor  in  the  denominator,  which  has  been  called  a  correlation  or 
coherency  factor  (3],  is  critically  dependent  on  the  phase  difference  between 
the  signals.  In  the  most  favorable  situation  *e(ye^)  ■  0  and  no  degradation 
is  caused  by  not  knowing  the  direction  of  the  interferer.  In  the  least 
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favorable  drcuaatance  Sc  “  |y|  and  tha  correlation  factor  la  Hy|2 


From  theae  axaaplea  w  conclude  that  the  efficiency  factor  la  due  to  the 
unknown  strength  of  the  Interference,  while  the  coherency  factor  la  due  to 
unknown  direction. 

When  Multiple  observations  are  available,  it  la  Instructive  to  compare 
the  reduced  Information  nat rices  for  the  Incoherent  and  coherent  cases.  For 
the  Incoherent  case,  we  have  from  (2.15) 

,<-)  .  R.  (Rj"  'j  [£,) 

(2.20) 


_  2r,(2wl)i 

- 7T- 


symmetrlc 


Re{Y  l  pnlpn2} 

r  M2 


The  bound  depends  on  the  actual  set  of  complex  signal  amplitudes  which 
occurred.  If  the  two  signals  are  really  Incoherent,  the  off-diagonal  term 
will  be  small  and  the  errors  will  be  almost  Independent. 

The  completely  coherent  case  gives 


F 


(-) 

uu 


2n(2vl)2  .2 
- - - d 

a 


i  i2  * 

IPjj  Retype) 

symmetric  fp2|2 


(2.21) 


The  coupling  tern  here  Is  not  reduced  by  multiple  observations,  and  so  the 

2  la 

resulting  bound  on  the  variance  contains  the  correlation  factor  1-Re  (yeJT). 


One  Emitter  In  the  Presence  of  One  Interfering  Emitter  and  White  Noise 

The  amplitude,  phase,  and  direction  of  the  emitter  are  modelled  as 
unknown  parameters.  The  Interference  Is  modelled  as  a  signal  from  a  known 


16 


direction  Uf  with  •  coaplex  Gaussian  amplitude  of  known  power.  The  covariance 
■atria  of  Interference  plus  noise  is  thus 

A  ■  g2(I  +  p  j  wt  Vj) 


with  Inverse 


where  y 


,-1  1  /T  H. 

A  ■  -j  (I  “  Y  v*  Y-i' 
a 
Pt 

— tt  .  After  considerable  algebra  we  find 

pj+l 

2  g2  1 _ 1 _ 

°u  |pn|2  (2wi>2 _ YljjJ2 

(a**)2  inrlvjVgl2 


(2.22) 


This  is  the  same  as  (2.18)  except  for  the  factor  y  >  which  approaches  1  as 
the  Interference  becoaes  large  relative  to  the  noise.  Note  also  that  as  the 
Interference  power  goes  to  sero,  (2.22)  becoaes  the  one  eaitter  bound 
(2.17).  This  bound  will  be  coapared  later  to  a  aiallar  one  in  which  the 
eaitter  signal  is  also  aodelled  as  a  Gaussian  process. 

F.  Craaer-Rao  Bounds  for  the  Gaussian  Signal  Model 

The  preceding  paragraphs  have  presented  Craaer-Rao  bounda  on  the  accuracy 
of  angle  of  arrival  estiaates  which  apply  when  the  signals  sre  deterministic 
with  unknown  parsasters.  Another  popular  signal  nodal  treats  the  coaplex 
amplitudes  as  coaplex  Gaussian  randoa  variables.  The  paraaeters  to  be 
estlaated  are  then  labedded  in  the  covariance  astrlx  of  the  observations 
rather  than  in  their  aean  value.  Zn  this  section,  CR  bounds  based  on  a 
Gaussian  signal  aodel  are  derived. 

Angle  estiaatlon  Is  to  be  bssed  on  H  Independent  snapshots  of  array 
data.  The  n*b  snapshot  la  of  the  fora 
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i.'J,  <W*V  +  £  .  •  Wl.H, 

Both  £  n  and  £  are  complex  circular  Gaussian  random  vectors,  l.e., 

"i  „>  •  '<1  „  I  •  0.  *<L 

8(£  D)  -  E<£  „  £  *)  -  0,  E<£  n  £  J)  -  P 


Consequently  z  Is  also  complex  circular  Gaussian. 


E  (z  )  -  E(z  i  V  0,  E(i  !  ")  *  VPV11  +  A  &  A 
—  n  —  n  —  n  —  n  —  n  0 


The  likelihood  function  for  the  observations  Is 
P<*  | .  •••«.„/«»?>•  "5T — 5T 

1  ”  v^lAl 


r  H  -1 
'I  *  n  A  i 
n 


» 

I 


with  logarithm 


In  p<£  ...  £  /  u_  ,  p)  ■  ~MN  Im  -  N  £n|A{  -  £  £ 

n 

Let  x,y  denote  any  two  of  the  unknown  parameters. 

A  typical  entry  In  the  Fisher  Information  matrix  is 


_  31 np  31np 
E  3x  3y 


The  derivative  is 


81  np 
8x 


n  LlAi_ 

|A  |  3x 


N 

l 

n*l 


H  3A 


-1 


3x 


Making  use  of  the  relations  [4] 
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m  e 


ft-  IMTrU 


3A  _  _A  -l  sa_  ,  -l 


this  becomes 


-  ? 

3x  n« 


M  H  3A  „  »  r.-l  3A, 

,llinirin-NTr{A 


-  N  Tr  A'1  A_1(S'A) 

1  N  H 

where  S  ■  jj-  \  £nAn  18  the  8anPJ-e  covariance  matrix. 

,  n**l 

we  have 


E  O  -  0. 


The  desired  expectation  Is 


(MaEjESE)  .  „  Tt  (A-1  fir1  |i) 

3x  3y  3x  3y 


Since  E(S)  -  A  , 


(2.23) 


This  result  Is  derived  In  detail  in  [5].  Appendix  2. 

Assume  first  that  the  "signal  in  space"  covariance  matrix  P  Is  completely 
known,  and  evaluate  the  information  matrix  for  estimating  £  .  The  required 
derivatives  are 


p  yH  3v_ 

3u±  3ut  3ua 


and  since 


3  v(u. ) 


|i-  -  (0  ,  ....  0,  -|-i-  .  0,  ..0)  ■  v(Ul)  t  .  J 


where  «<  li  a  unit  vector  whose  1th  element  Is  unity,  and 


. -ar* 
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we  have 


V  -  ..  £(11^)) 


3A  •  H  „  „H  .  H  •  H 

3u^  i  —  1.PV  +VPil£lV 


The  ik^  element  of  the  information  matrix  ia  given  by  (2.23)  aa 


Flk  -  E  (M2E  MSR  , 


3u_,  3^ 


N  Tr{A_1(V  e  ±  e  ®  PVH  +  VPe  ±  e  ^VH)  A_1  e  fc  e  ®  PVH+  VPe  fc  e  ^H)> 


As  in  Section  C,  we  define  the  matrices 
W  -  V11  A-1  V 

w  -  vha-1v 

W  -  V^A_1V 

Since  any  cyclic  permutation  of  matrlce8  does  not  affect  the  trace  operatic#, 
the  preceding  expression  can  be  rewritten  aa 

Fik  "  N  Tr£(®  ^  PW  e  k>  (e^  k  PW_ei)  +  complex  conjugate 
+(<e  ^  W  e  k>  (£  k  PVfPe  i>  +  complex  conjugate) 

The  quantltiea  in  parenthe8la  are  scalars,  so  the  trace  operator  can  be 
j  dropped.  Also  each  term  can  be  combined  with  its  complex  conjugate  to  give 

Fik  -  2  N  Re  {(e  ”  PW  e  k)  (e  J  PW  e  J  +  (e  J  W  e  fe)  (e  J  PWP  e  t>} 

j 


20 


Define  the  Hadaaard  product  of  two  MXN  aatrlces  to  be 


(A  O  B)  -  A  B 

win  Bn  i 


Then  we  can  write  the  natrix  [F^,]  as 


P  -  2N  Re  (P&  <0  (P&)T+  W  O  (PWP)T> 


u  u 


The  result  can  be  put  in  a  wore  convenient  fora  [3]. 
A  of  the  fora 


(2.24) 


Seek  an  inverse  for 


r1  -  A-1  -  A-1  V  Q  V®A 


-1 


Coaputatlon  yields 

AA_1  -  I  +  V(P-Q-PW  Q) V^a""1  -  I 
o  o 

A_1A  -  1  +  a“1V(P-Q-QM^P)VH  -  I 


where  W  -  V^A^V. 
o  o 


Thus  Q  must  satisfy 


P-Q-PV  Q-QW  P 
o  o 


Solving  for  Q,  we  find 


Q  -  (I  +  PW  )_1  P-P(I  +  W  P)'1 
o  o 


-1 


and,  if  P  exists 


q  -  (p"1  +  Wo)_1 


(2.241) 
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Computation  of  the  elements  of  (2.24)  with  frequent  use  of  (2.241)  yields 


F  ■  2N  Re  {H  -  )  O  (P-Q)T  +  qft  O  (Q&  )*>  (2.242) 

o  o  o  o  o 


where 


&  -  v^“H 

o  o 

H  -  1  V 

o  o 

This  Is  the  desired  result  for  the  case  of  Gaussian  signals  with  known 
covariance. 

As  the  smallest  (diagonal)  element  of  P  becomes  large.  Q  ♦  ,  P-Q*-P 

and  (2.242)  becomes  asymptotically. 

F  -  2N  Re((W  -  ft8  H_1  &  )  O  P)  (2.243) 

o  o  o  o 

In  the  case  of  unknown  covariance  P,  Indexing  becomes  a  problem.  He 
shall  restrict  our  attention  to  one  simple  limiting  case,  namely,  completely 
incoherent  signals  (diagonal  P).  He  again  make  use  of  (2.23).  The  necessary 
derivatives  are 


wr m  dpf:  {Ao + 1  (u»^H(u«)} 

11  11  m 


IT  H 

VS  it  t  v 


He  then  compute 


or 


E(H2E  >  *  N  Tr{A_I(V  e  e  ?  VH)  A-1(V  e  .e  ?V8)} 
3Pli  3Pkk  ~  k~  k 

-  N  |e  J  VHA_1Ve  J2 

E(2||E)Ti||E.N(W0HT) 
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4S*  :  ?v- 


WH2*  -  N  Tr{A_1Ve  fe  JvV^Ve  ^  J  PVH  +  VPe  k«  JV0} 


2N  Ref (e  JvVHe  fc)  (e  [W^Ve  ±)} 


or 


B{  -  2N  Re  {&  O  (PW)T) 


The  overall  lnforaatlon  matrix  for  a  diagonal  covariance  matrix  P  then  takes 
the  form 


F  F 

P  P  P  u 

w  o  WT 

2  Ref 6  O  (PW)T) 

F  - 

F  F 

U  P 

-  N 

symmetric 

2  RefW  O  (PWP)T  ♦  pfr  O  (PW)T) 

m  m 

• 

. 

(2.25) 


Uae  of  (2.241)  gives  the  more  convenient  form 


F-N 


0 

I 


QW  O  (QW  )T 
o  o 

symmetric 


2  Re{Q&  O  (QW)1} 
o  o 

2  Ref(V  -W®q6  )  O  (P-Q)T+QW  O  (QW  )T} 
o  o  o  o  o 


o 


o 

i 


(2.251) 


provided  that  P-*  exists. 


G.  Examples  of  CR  Bounds  for  Gaussian  Signal  Model 
One  Emitter  in  White  Noise 


The  covariance  matrix  for  this  case  la 


2  n 

A  ■  o  (I  +  p  v  v  ) 


where  p  la  the  array  slgnal-to-nolae  ratio  per  obaervatloo.  The  lnveree  la 


.-1  1 


(I - |j-r  v  vH) 

P+1 - 


We  compute 


W  -  vHA  _1v  ■  —  (1 - £— )  •  - - - 

-  -  2  p+l'  1.  . . 

00  (p+1) 

t.  H.-l»  1  ,,  p  x  H*  1  J  r  *■ 

W.vA  V  ■  — 2  (1  -  v  V  -  -5 - i  I  j— 

o2  p  1  a2  (p+1)  M  X 

V.  »H  .  -1  •  1  •  H*  p  i  H»,2 

W-vA  V  •  -y  V  V  -  -gr* -  |  v  vj 

O  O  (p+1) 

We  assume  as  before  that  the  array  coordinate  system  Is  chosen  so  that 

l  *.-»  • 


r  (^.(W 


Equation  (2.2A)  becomes 


F  -  2N 
uu 


(op)2  (2*1)2  _  2ty  ,2x1x2 

_2/_,x  _2  I+l7p  1  X  ' 


o  (p+1)  o 


Therefore  the  CR  bound  for  estimating  the  array  coordinate  u  of  one  Gaussian 
signal  form  N  independent  observations  is 


,2  >_!+!/£_ 

U  2Np(2*tr 


(2.27) 


which  is  the  same  as  the  bound  for  a  deterministic  signal  (2.17)  except  for 
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This  differs  froa  «  previous  example  In  that  the  desired  signal  Is 
node lied  as  Gaussian  rather  than  deterministic. 


The  covariance  matrix  is 
A  -  R  +  P  v  vH 

where  R  Is  the  known  covariance  aatrix  of  the  Interference  and  nolee  and  P  is 
the  unknown  power  of  the  emitter.  For  brevity,  define 


H  -1 

r  ■  v  R  v 


•  H  -1  • 
r  ■  v  R  v 


•1L-1  • 
r  *  v  k  v 


Then 


The  required  matrices  (scalars)  are 


W 


r 

1  +  Pr 


r 

I  +  Pr 


W  -  r 


P 

1  +  Pr 


|r| 


2 


Using  (2.25)  we  obtain  the  Information  matrix  for  the  paraaetere  P,  u 


respectively. 


„  f  r2  2Pr  Re(r)  1 

_w  .  |(2.5 

(1+Pr)2  symmetric  2P2((H-Pr)(rr-|r|2>  ♦  2  Re2(r)]| 


This  agrees  vi-n  a  result  obtained  by  Miller  and  Huber  (6) • 
The  reduced  information  matrix  for  u  alone  la 


•r  -  *  &  <-  -  w*) 


(2.29) 


To  compare  this  with  a  previous  result  (2.22)  which  was  derived  using  a 
deterministic  signal  model,  we  take 

R  -  a2  (I  +  px  Vj  vj> 

The  Interference  then  consists  of  white  noise  plus  a  single  source  of  GauBslan 
Interference  at  known  direction  Uj.  Then 

R~l  -  \  (X  -  yvjvj) 

o 

where  y  *  Pj/tpj+l). 

We  evaluate 

H  -1  1  ,,  ,  H  ,2. 

r  ■  ^  R  -  —  (1-YlVjvJ  ) 

a 

r  *  i.  R  \ - 2  (Vi>  '-V.* 

0 

; .  ttjli  -  -l 

a  a 

Substitution  of  these  results  Into  (2.29)  gives 


(2.30) 


2p#[(2irl)2 


1  + 


P.Ci-y  |v*vj2 


] 


p 

where  PB  “  “J  • 

<j 

The  reciprocal  of  (2.30)  agrees  with  (2.22)  except  for  a  factor  (the 
denoalnator  of  (2.30))  which  approaches  unity  for  large  signal-to-noise 
ratio. 


i 
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2.  ZIV-ZAXAI  BOUNDS  1 7 ] 


A.  One  Unknown  Parameter 

Let  u  be  an  unknown  pereaeter  taking  values  in  an  interval  U.  The  value 
of  u  is  to  be  estiaated  baaed  on  soae  set  of  observations  which  aay  be  written 
as  a  vector  t,  A  statistical  model  is  given  relating  the  observations  to  the 
value  of  u. 

For  any  two  values  of  u,  uj ,  and  U2  >  uj ,  there  exists  an  optlaal 
statistical  decision  rule,  the  likelihood  ratio  test  (LRT),  for  deciding  on 
the  basis  of  the  observations  which  of  the  two  values  is  correct.  The  LRT  is 
optlaal  in  the  sense  that  it  minimizes  the  probability  of  error. 

If  the  values  Uj ,  u2  are  equally  likely  to  occur  a  priori*,  the  LRT 
chooses  uj  if 


,  e'ilV 
'  -‘“HIIsT 


and  otherwise  u2.  The  probability  of  error  for  this  optimal  test  is 

Pg  -  yPr  (i  <  0  |  uL)  +  |Pr  (1  >  0  |  u2) 

Note  that  PE  depends  on  the  two  values  Uj,  u2,  and  also  on  the  statistics  of 
the  observations. 

Now  consider  a  suboptimal  decision  rule  based  on  an  arbitrary  estimate  u 

A 

of  u.  This  rule  chooses  u^  if  u  is  closer  to  uj  than  to  u2,  and  vice  versa. 
The  error  probability  of  this  decision  rule  is 


It  is  possible  to  derive  a  more  general  fora  of  the  ZZ  bound  by  allowing 
arbitrary  a  priori  probabilities  [8].  However,  this  tends  to  obscure  the 
basic  ideas  Involved  and  furthermore  provides  no  iaproveaent  in  the  bound  for 
the  examples  to  be  discussed  here.  See  Appendices  II,  III. 


Application  of  the  Chebychev  Inequality  ylelda 


e2(u1> 


u,  -  u. 

(-V) 


-_  +  i. 
2  2 


e2(u2) 


u«  -  u, 

(-V1) 


2  * 

where  e*(u)  denotea  the  mean  aquared  error  in  the  estiaate  u  when  u  ia  the 

true  value.  Since  Qg  >  Pg,  we  have 


J  (e2(uj)  +  e2(u2))  >  (-2-^ — -)  Pg(uj ,  u2> 


(2.31) 


for  any  pair  of  valuea  Uj,  u2. 

The  right-hand  aide  of  (2.31)  haa  been  increaaed  by  a  factor  of  2  by  Vax 
and  Ziv  [9]  uaing  a  such  wore  complicated  argument  (aee  Appendix  I).  We 
Incorporate  thla  improvement  in  the  remainder  of  thia  report.  Thua ,  the  baaic 
Ziv-Zakal  bound  ia 

2 

19  9  ^  ^9  ^  |  ) 

y  (•  (ut)  +  *(u2))  >  -  ■  y  -  PE  (ur  u2)  (2.32) 

Observe  that  the  ZZ  bound  la  a  lower  U>und  on  the  average  mean  aquared 
estimation  error  at  two  parameter  values.  Thla  is  the  beat  one  can  hope  for 
when  considering  a  completely  unrestricted  class  of  estimators,  since  the 

A 

estimate  u  ■  u^  haa  aero  error  when  u  ■  uQ.  To  obtain  a  point  bound,  note 
that  for  any  given  estimator  u  there  exists  some  value  u«  which  produces  the 
largest  mean  square  error.  The  left  hand  side  of  (2.32)  la  then  less  than 
e2(u*)  for  any  choice  of  uj,  u2  and  so  we  can  maxi  mice  the  right-hand  side 
over  these  variables,  getting 

29 


(2.33) 


2 

e  (uj  >  max 

*  V  UJ 


(u  -  u  Y 

-■  2  —  PE  <U1*  U2) 


which  la  a  lower  bound  on  Che  worat  poaaible  error  any  eatimator  can  make. 


B.  Additional  Unknown  Parametera 

The  argument  leading  up  to  the  two-point  bound  (2.32)  la  unchanged  by  the 
presence  of  other  unknown  parameters  In  the  statistical  model*.  Let  a,  denote 
the  vector  of  additional  parameters.  The  LRT  chooses  between  two  values  (u^ , 
oj),  (u2,  a^)*  I ts  error  probability  la  a  function  of  these  values. 

As  before,  define  a  suboptimal  hypothesis  test  based  on  an  arbitrary 

A 

estimate  u  of  u.  This  leads  directly  to  a  multi-parameter  version  of  (2.32). 


I  .  2  (u9  ^ , ) 

j  (e  (uj,  Oj)  +  e  (u2,a2))  >  - g - PE^U1’  u2*  —1*  —2^  (2*3^) 


There  are  many  ways  to  obtain  point  bounds  from  (2.34).  The  parameter  vector 
a.  can  be  divided  into  a  vector  of  "interesting"  parameters  and  a  vector  of 
"nuisance”  parameters  y_.  The  distinction  is  that  we  wish  to  obtain  a  bound 
which  is  an  explicit  function  of  the  "interesting”  parameters,  but  are  willing 
to  choose  worst-case  values  y^Cfl)  for  the  "nuisance"  parameters  (and  u).  This 
gives 


(6.J  ))  +  e^u*(B.2)42,x*(i2»  J 


(u2~u1) 

>  _  P£ (u  J  ,u2  ,J6  J  ,6.2  »12 ) 


Known  parametera  are  simply  part  of  the  statistical  model  for  the 
observations  and  need  not  be  considered  explicitly. 


The  inequality  is  valid  for  any  choice  of  Jj,  Qj ;  however,  it  is  most 
useful  when  we  choose  dj  “  j*2  *  £  to  get  a  point  bound 


2  <u»  -  u.)2 

el  (u*(B),  £,  J[*(8))  > - 5 - PE^U1*  u2*  ll»  I** 

u 

This  inequality  holds  for  any  choice  of  uj,  U2,  yj,  yj ;  as  in  the 
one-parameter  case,  the  tightest  bound  is  obtained  by  maximizing  the  right- 
hand  side  with  respect  to  these  variables.  The  result  is 


e2(u*(6),  6,  x*(£)) 
u 


> 


l“i2  Pe<ui'  V  *!• 

(2.35) 


The  problems  to  be  considered  in  this  report  are  signal  parameter  esti¬ 
mation  problems  characterized  by  the  observation  model 


K  j4»t 

z  -  l  A  e  v(u. )  +  £ 
i-1  1 


To  get  a  meaningful  bound  on  the  error  in  estimating  u^,  the  associated 
signal  amplitude  A*  must  be  treated  as  an  "interesting"  parameter.  The 
reason  Is  that  the  worst-case  value  of  A^  is  clearly  zero.  The  signal  is 
then  absent  and  the  observation  provides  no  information  about  u^.  The 
other  amplitudes,  however,  may  be  treated  as  "nuisance”  parameters. 

The  Cramer-Rao  bound  clearly  shows  how  the  presence  of  additional  un¬ 
known  parameters  degrades  the  bound  on  accuracy  for  estimating  a  parameter. 
With  the  Zlv-Zakal  bound,  the  distinction  is  less  clear.  The  parameters  £  in 
(2.35)  are  held  fixed  in  the  two  hypotheses  being  tested;  they  are  thus 
treated  in  the  same  way  as  the  known  parameters.  The  parameters  are 
different  in  the  two  hypotheses  being  tested  and  therefore  wst  be  unknown 


parameters.  The  bound  (the  right-hand  aide  of  (2.35))  is  an  explicit  func¬ 
tion  of  parameters  £  and  makes  no  distinction  between  whether  they  are  Itnown 
or  unknown;  it  is  an  implicit  function  of  unknown  parameters  y_  for  which 
worst-case  values  have  been  selected. 

We  now  specialize  the  bound  (2.35)  to  some  cases  of  particular  inter¬ 
est.  Paralleling  the  development  of  Cramer-Rao  bounds,  we  begin  with  deter¬ 
ministic  signal  models  and  then  treat  the  Gaussian  signal  model. 

C.  One  Known  Signal  in  White  Gaussian  Noise 

The  observation  is  a  complex  Gaussian  random  vector  with  covariance  ma¬ 
trix  o^I  under  either  hypothesis;  its  mean  is 

Ae^jKu)  . 

The  error  probability  for  deciding  between  two  values  u}  and  U£  when 
A  and  $  are  known  is  (see  Appendix  II)1 . 

PE  -  erf  c^  [  —  |  s(u^)  -  s(u2)  j] 

Of  particular  interest  is  the  case  where  the  length  of  the  signal  vector  £  is 
unity  for  any  parameter  value,  and  the  projection  sfu^^sf^)  depends 
only  on  the  parameter  difference  u^  -  U2»  These  assumptions  hold,  for 
example,  in  the  tlme-of-arrival  and  direction-finding  problems  when  only  one 
signal  is  present. 

Then  the  error  probability  becomes 


PE  ’  er£c*  t  £  /  1  -  Re{c(uj  "  u2)K 


erfc*(x)  ^  --  1 -  /  e  2  dr 

/  2v  a  x 
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where  |js(u)|2  ■  1  and  c(u^  -  u, )  ■  £(Uj)Ws(u2)«  The  function  c  is  a  measure 
of  the  correlation,  or  degree  of  similarity,  between  normalized  signals  hav¬ 
ing  different  parameter  values.  The  ZZ  bound  for  this  case  is  then 


V  [e2(ulfA)  +  e2(u2,A)]  2 


(U1  "  u2> 


irfcJ« 


1  -  Re{ c(u  ^  -  u2>] 


(2.36) 


and  the  worst-case  Z  bound  is 


e_(u4.,A)  >  max 

~  U1  '  u2 


(U1  "  U2r 


erfc*(  —  /  1  -  Re{  c(u,  -  u» )] 


(2.37) 


The  right-hand  side  of  (2.37)  involves  two  competing  terms  -  the  first  grows 
quadratically ,  while  the  second  usually  becomes  small  as  the  parameter  values 
separate.  In  the  time-of -arrival  problem,  for  example,  c(uj  -  U2)  is  the 
normalized  correlation  function  of  the  pulse  envelope,  which  typically  de¬ 
creases  smoothly  and  becomes  zero  once  the  time  difference  uj  -  U2  ex¬ 
ceeds  the  pulse  width.  Typical  behavior  of  the  right-hand  side  of  (2.o7)  is 
shown  as  a  function  of  parameter  separation  for  increasing  values  of  signal- 
to-noise  ratio  in  Fig.  2.1.  The  non-monotonic  behavior  reflects  that  of  the 
signal  correlation  function  c. 

At  large  signal-to-noise  ratios,  the  maximum  required  on  the  right-hand 
side  of  (2.37)  occurs  when  uj  -  U2  is  small,  as  shown  in  Fig.  2.2.  The 
argument  of  the  error  function  can  be  bounded  for  small  e  *  uj  -  U2  by 
using  the  inequality  cos  x  ;•  1  -  0.5  x-. 


£  V1  -  Re-i  '  e 

0  »  M  - 


j2*t-  e 


2itx  c 


!  -f V1  C08  2"£r ^V1  -  s*  r -)  ) 


2r  l  t 
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Behavior  of  Ziv-Zakai  bound  vs.  emitter  separation  for  increasing 


f  Zlv-Zakai  bound  for  large  S/N  ratio. 


where 
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The  derivation  of  the  ninitaun  error  probability  for  deciding  between 
Hi  and  H2  is  given  in  Appendix  III.  Assuaing  |s(uj.)|2  ■  1, 

£h(ui)  £(u2>  ■  c(ui  -  U2>  as  in  the  preceding  section,  the  result 

is 

PE(ul’  u2}  "  erfc*^/  1  -  | c(Uj  -  u2)|j 
and  the  corresponding  worst-case  bound  is 

e2(uA)  2  — ~~2 — ” —  erfc*(£  /  1  -  |c(u  -  u2)|)  (2.39) 

~VU2 

Thus,  the  only  change  to  the  bound  (2.37)  is  that  Re  c  is  replaced  by  |cj. 

Note  that  by  treating  phase  as  an  unknown  (nuisance)  parameter,  we  ob¬ 
tain  a  bound  which  differs  from  that  for  known  phase  and  is  independent  of 
the  actual  phase  value. 

It  was  pointed  out  previously  that  at  large  S/N  ratios,  the  desired 
maximum  occurs  for  snail  Uj-Uj.  this  laPlles  |c(u^-u2)|  '  Re  c(Uj-u2)  ”1. 
Thus,  the  worst-case  ZZ  bounds  for  known  and  unknown  phase  coincide  for  large 
S/N  ratios. 

E.  One  Signal  with  Random  Phase  in  White  Gaussian  Noise 
A  third  signal  model  of  interest  is  one  that  treats  phase  as  a  random 
variable  uniformly  distributed  as  (0,  2x).  It  can  be  shown  [8]  that  the  op¬ 
timum  hypothesis  test  between  equally  likely  parameter  values  u}  and  U2 
has  error  probability 


P  -  S(a,B)  -  |  (l  -  Q  (/0,  ✓£))  +  \  Q(/o,  /0  ) 


where 


2a 


(  1  -  V  1  *  c<ui  "  u2>|  ) 


6  -  •—  (  i  +  V1  "  |c(ui  _  u2} | 2) 

2a 


and 


2  2 

x  +  a 


Q(a,  b)  -  /  xe 
b 


IQ(ax)  dx 


Is  Marcum's  Q  function.  The  resulting  worst-case  ZZ  bound  is 


2  ^U1  u2^ 
e  (u*)  >  max  - ^ -  S(o,  0) 

ul““2 


(2.40) 


For  large  S/N  ratios,  the  expression  for  the  error  probability  is  asymptoti¬ 
cally  equal  to  erfc*  (/$  -  /oT  )  and  the  right-hand  side  of  (2.40)  becomes 
(approximately ) 


max 

u 


r erf V 1  +VT“-  |c(u)"p  -  y i-  V  1  *  |c(u)|i 
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At  high  S/N  ratios,  the  maximum  occurs  for  small  u,  so  that  |c(u)|  -  1,  and 
the  above  expression  is  approximately 


Thus,  at  high  S/N  ratio,  the  random  phase  bound  agrees  with  the  unknown  phase 
bound. 

F.  Multiple  Elementary  Signals  with  Unknown  Phases  in  White 
Gaussian  Noise 

The  observation  is  modelled  as 
£  -  V(u)  £  +  £ 

The  columns  of  V  are  the  elementary  signals  v(u^)  and  £  is  a  vector  of 
unknown  complex  amplitudes.  The  minimum  error  probability  for  deciding 
between  (u,£)  and  [u* ,  £'  ]  in  this  case  is  derived  in  Appendix  IV.  The 
result  is 

PE  "  ) 

✓  2 

where 

d  -  f(u,  u’ ,  £,  £' ) 

The  resulting  worst~case  ZZ  bound  on  the  error  in  estimating  the 
angle-of -arrival  un  of  the  nth  signal  is 

-  (u  -  u’)2  Ad 

e  (u  A  )  >  max  - 5 -  erfc*(-£-  )  (2.41) 

”  V2V 

The  maximization  over  p,  p'  can  be  done  analytically,  as  discussed  in 
Appendix  IV.  However,  the  maximization  over  u,  u'  does  not  appear  tractable 
for  more  than  one  signal. 


G.  One  Gaussian  Signal  In  White  Gaussian  Noise 


The  contending  hypotheses  are 


»x  !  i„  ■  *<“l>  + 

®2  !  ‘  'n  2<02>  +  i. 


n  ■  1,  ...N 


where  ^  is  a  complex  Gaussian  noise  vector  with  zero  mean  and  covariance 
o 

matrix  a  I,  and  pQ  la  a  complex  Gaussian  random  variable  with  zero  mean  and 
variance  MP.  Both  ^  and  pn  are  independent  from  look  to  look. 

The  minimum  error  probability  for  deciding  whether  or  H2  Is  true  Is 
shown  In  Appendix  V  to  be 


P£  -  Pr  {U  >  V} 


where  U,  V  are  Independent  chl-squared  random  variables  with  2N  degrees  of 
freedom  and  variances  (per  degree  of  freedom) 


2  *  H  1  +  P I vf  v- | 2  , 

a v  ”  P  (1  +  2Re(A  v®^)  + - r"+"p -  ial  ^ 

2  ,i  ,2  .  *n  ,  .  1  +  p  1-1-2^ 


p  (  | a |  +  2Re(A  v^)  + - 


) 


and 


P 

a 


X 


MP 


~2  m  array  slgnal-to-nolse  ratio/look 


X 


-  AhTT 


2  44 


.  ,  ,  H  ,2 

1  + 


1  + 


1  +  P 


It  is  further  shown  In  Appendix  VI  that  this  probability  is  equal  to  the 
probability  of  obtaining  less  than  N  successes  in  2N-1  trials  with  probability 


of  success 


P 

8 


2.2 
a  +  o 
u  v 


The  resulting  worst  case  ZZ  bound  Is 

2  <ui~“2>2 

•  <u)  *  Sjs»2  "i  —  PE 

It  Is  particularly  Interesting  to  look  at  the  result  for  a  single  snapshot 
(N-l).  Then 


A  rather  tedious  calculation  shows  that  for  large  slgnal-to-nolse  ratio  p 
Pg  •  [(1  -  | c<u) | 2 )  pi  1 


and  the  ZZ  bound  becoaes 


2,  .  „  1  u 

e  Cu)  >  max  — 


u 


P  i  - 


ic(«>r 


This  Is  quite  different  from  the  corresponding  results  for  deterministic 
signal  models.  In  those  cases,  the  maximizing  value  of  u  exhibits 
discoutlnultles  as  slgnal-to-nolse  ratio  Is  increased,  and  ultimately 
approaches  zero.  In  this  case,  the  maximizing  value  of  u  is  independent  of 
slgnal-to-nolse  ratio  (provided  S/N  Is  large  enough)  and  is  not  necessarily 
small.  Thus  the  ZZ  bound  does  not  approach  the  CR  bound  at  high  slgnal-to- 
nolse. 


As  the  number  of  looks  Increases,  the  behavior  of  the  Gaussian  signal  ZZ 


bound  approaches  that  of  Its  deterministic  counterpart.  Fig.  2.3  shows  the 
two  sets  of  bounds  for  a  uniform,  9-element  linear  array  with  X/2  spacing,  for 
1,  5,  10,  and  100  looks.  The  two  bounds  ultimately  agree  for  2  or  more  looks; 
with  20  looks  they  are  almost  Identical  for  all  signal- to- noise  ratios. 

The  Intuitive  explanation  of  this  result  Is  that  on  a  single  trial  the 
probability  of  drawing  a  low  slgnal-to-nolse  ratio  from  a  Gaussian  sample  Is 
too  high.  Ulth  multiple  Independent  observations,  this  situation  Is 
rectified. 

The  conclusion  Is  that  when  multiple  Independent  observations  are 
available,  the  choice  of  a  signal  model  is  not  critical.  However,  If  only  one 
observation  Is  available,  the  two  signal  models  lead  to  quite  different 
results. 


ARRAY  S/N  RATIO  PER  LOOK  HB) 


III.  UNEQUALLY  SPACED  LINEAR  ARRAYS 


A.  Applicable  Methods 

Many  angle  estimation  techniques  can  be  interpreted  in  terms  of  adaptive 

array  concepts.  The  basic  idea  is  to  form  a  weighted  sum  of  the  array  element 

outputs  so  that  the  total  output  power  is  minimized.  Some  sort  of  constraint 
must  be  imposed  on  the  weighting  vector  to  avoid  the  trivial  solution  v_  ■  0. 

If  the  constraint  is  a  linear  one  of  the  form  -  1.  then  it  is  easily 

shown  that  the  minimizing  weighting  vector  is 


and  the  adapted  array  pattern  in  the  direction  u  ■  sln8  is 
G(u)  •  |wHv(u)J  *  (c^R  \v(u)| 

When  the  signal-to-noise  ratio  (SNR)  is  large  enough,  it  can  be  shown  that 
under  certain  conditions  the  adapted  pattern  has  (asymptotic)  nulls  in  the 
directions  of  the  external  sources.  The  required  conditions  are: 

(1)  The  direction  vectors  of  the  sources  are  linearly 
independent . 

(2)  The  signal-in-space  covariance  matrix  P  has  an  inverse. 

These  nulls  can  be  used  to  determine  the  source  directions. 

If  £  ■  (I,  0,  ...,  0),  the  constraint  fixes  the  weight  applied  to  one 
element  (the  "reference"  element).  For  a  uniformly  spaced  linear  array,  this 
method  is  mathematically  Identical  to  a  special  case  of 
autoregressive/maximum-entropy  spectral  estimation  in  which  the  model  order  is 
equal  to  the  number  of  array  elements.  While  it  continues  to  have  meaning  in 
an  adaptive  array  context  for  unequally  spaced  arrays,  it  does  not  provide  a 
spectral  estimate  having  maximum  entropy.  Since  the  weights  constitute  the 
best  linear  prediction  of  the  signal  at  the  reference  element  based  on  the 
signals  observed  at  all  the  other  elements,  we  will  refer  to  this  method  as 


the  linear  prediction  (LP)  technique.  In  adaptive  array  terminology,  when  the 
reference  element  is  the  main  antenna  and  the  other  samples  are  from 
auxiliary  antennas,  the  technique  is  called  sidelobe  cancellation. 

When  £_  ■  _v(uq)  ,  the  constraint  fixes  the  adapted  array  gain  (w^v(u)|  in 
the  direction  uq.  This  pattern  has  (asymptotic)  nulls  In  all  source 
directions  (except  ug),  which  could  be  used  to  locate  them.  However,  the 
rationale  behind  the  maximum  likelihood  method  (MLM)  is  somewhat  different. 

Assume  that  the  process  on  the  aperture  consists  of  a  desired  signal  from 
direction  uq  and  interference.  Assume  further  that  the  signal  and 
interference  are  uncorrelated.  The  expected  power  output  of  the  array  using 
weighting  vector  Is 

E  |wHzJ2  -  A  w 

where  R  -  RT  +  P  v(u  )  v  (u  )  is  the  covariance  matrix  of  the  process.  The 
constraint  w_v(u0)  “  1  fixes  the  output  power  due  to  the  desired  signal. 
Minimizing  the  interference  power  output  is  thus  equivalent  to  minimizing  the 
total  power  output.  The  minimizing  weighting  vector  is 


R_1v(un) 

w(u  )  ■  ■■  *  ■  ■ 

-o'  H.  S  D-1  ,  \ 

v  (u  )  R  v(u  ) 
—  o  —  0 


and  the  expected  output  power  is 


p(uo>  -  “g - *=1 - 

vH(u  )R  v(u  ) 
—  o  —  o 


(3.1) 


This  quantity  is  computed  as  a  function  of  uQ  and  the  locations  of  Its  peaks 
(exceeding  a  threshold)  are  the  estimates  of  the  source  directions. 

A  third  choice  for  the  weight  constraint  Is  ■  1.  This  fixes  the 
thermal  noise  output  power  of  the  array.  The  minimizing  weighting  vector  In 
this  case  satisfies  the  equation 
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Rw  «■  Xw 


r 


This  eigenvalue  problem  has  M  solutions,  M  being  the  number  of  array  elements 
and  consequently  the  dimension  of  R.  The  M  solution  vectors  are  mutually 
orthogonal  and  their  associated  eigenvalues  are  real.  The  desired  solution 
is  the  eigenvector  corresponding  to  the  minimum  eigenvalue.  In  general,  this 
solution  is  not  unique.  It  can  be  shown  [101  that  if  I  <  M  emitters  are 
present,  iR  has  I  "signal”  eigenvectors  with  relatively  large  eigenvalues  and 
M-I  "nois4*"  eigenvectors  with  a  common  minimum  eigenvalue.  Any  one  of  these 
1  noise  eigenvectors  is  a  valid  solution  to  the  minimization  problem. 

If  the  signal-in-space  covariance  matrix  P  ia  non-singular,  it  can  be 

shown  that  the  signal  vectorib  v(u^),  i  "  1,  •  •  •  I  and  the  "signal" 

eigenvectors  span  the  same  I-dimensional  subspace,  which  is  orthogonal  to  the 

(M-I)-dimensional  subspace  spanned  by  the  "noise"  eigenvectors  It  follows 

H 

that  the  adapted  array  patterns  w^  v(u)  all  have  nulls  in  each  of  the  signal 
directions.  When  the  covariance  matrix  is  perfectly  known,  these  are  true, 
rather  than  asymptotic,  nulls. 

The  array  pattern  generated  by  a  particular  "noise"  eigenvector  will,  in 
general,  have  extraneous  nulls.  However,  the  only  nulls  common  to  all  of  the 
■  "noise"  eigenvector  patterns  are  those  corresponding  to  source  directions. 

This  is  true  because  if  _v(u)  is  orthogonal  to  all  the  "noise"  eigenvectors,  it 

must  lie  in  the  signal  subspace,  which  is  spanned  by  the  source  direction 
i  vectors  jKuj),  .  .  .,  _v(uj).  But  v(u)  is  linearly  independent  of  these 

vectors  unless  u  Is  one  of  the  source  directions  .  The  appropriate  function 

11  to  examine  for  nulls  is  thus 

) 

i 

! 

♦ 

: 

*If  the  direction  vectors  are  not  linearly  independent,  extraneous  nulls  will 
|  exist. 
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for  which  all  nulls  correspond  Co  source  directions.  This  technique  Is  called 
MUSIC  (Multiple  signal  classification)  [10). 

As  the  signal-to-nolse  ratio*  becomes  large,  it  was  shown  in  [3]  that 


where 


EN  -  w2“-**  -Sm-iI 

is  a  matrix  whose  columns  are  noise  eigenvectors.  The  MLM  method  seeks  the 
maxima  of  (3.1),  or  equivalently  the  minima  of 

vH(u)  R  1  v(u)  ♦  |E^  v(u)|  -  l  |wjj|  _v(u)|2 

Thus,  MLM  is  asymptotically  (high  S/M)  equivalent  to  MUSIC  and  has  no 
extraneous  nulls.  At  lower  S/N,  however,  extraneous  nulls  may  be 
encountered.  TMs  Is  also  true  for  the  linear  prediction  method. 

To  reject  whatever  extraneous  nulls  exist,  an  estimate  of  the  power 
received  from  each  candidate  source  direction  Is  made.  Those  that  fall  below 
a  threshold  are  eliminated.  The  power  estimate  is  obtained  by  solving  the 
matrix  equation 


R  *  V  P  +<j2I 

for  the  signal-in-space  covariance  matrix  P,  with  the  result 
P  -  V+  (R  -  ,2  I)  V+" 


Defined  as  the  ratio  of  the  smallest  non-zero  eigenvalue  of  P  to  the  noise 
power  o. 
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+  H  H  t 
where  V  ■  (V  V)  v  is  the  pseudo-inverse  of  the  matrix  of  (candidate) 

direction  vectors  V.  An  estimate  of  the  noise  power  level  is  required. 

The  rational  for  the  adaptive  methods  collapses  when  the  signal  sources 
are  coherent,  i.e.,  when  the  signal-in-space  covariance  matrix  P  is 
singular.  The  adapted  array  pattern  no  longer  has  nulls  in  the  directions  of 
all  sources.  The  reason  is  that  the  component  of  array  output  power  due  to  a 
combination  of  coherent  sources  can  be  reduced  to  zero  by  cancellation,  i.e., 
by  combining  them  with  appropriate  complex  weights.  These  weights  are  the 
values  of  the  adapted  pattern  at  the  directions  of  the  coherent  sources. 

A  technique  known  as  "spatial  averaging"  [3)  has  been  shown  to  be 
effective  in  locating  coherent  sources  with  a  uniform  linear  array.  At 
present,  no  such  method  is  known  for  unequally  spaced  linear  arrays. 

Extension  of  other  high  resolution  spectral  estimation  techniques  such  as 
maximum  entropy,  AR,  and  ARMA  modeling,  etc.,  to  unequally  spaced  linear 
and/or  two-dimensional  arrays  is  the  subject  of  much  recent  research  Ill, 
12].  We  plan  to  study  and  evaluate  some  of  these  algorithms  during  the  coming 
year. 


B.  Array  Geometry 

The  use  of  unequally  spaced,  small  linear  arrays  is  motivated  primarily 
by  the  desire  to  avoid  the  grating  lobes  which  result  from  uniform  element 
spacing*.  Elements  spaced  uniformly  d  wavelengths  apart  can  determine  the 
direction  u  *  sin6  of  a  source  only  modulo  This  results  from  the  fact  that 
phase  can  only  be  determined  modulo  2x  and  consequently  sources  at  any  two  of 
the  directions  un  *  sin60  *  nj  produce  identical  amplitude  and  phase  patterns 
at  the  array  elements.  This  is  no  longer  true  when  the  spacing  is  irregular. 


*Here  again,  we  assume  that  the  candidate  direction  vectors  are  linearly 
independent  (i.e,  the  direction  matrix  V  has  full  rank). 

*Unequal  spacing  is  also  used  as  a  form  of  tapering  in  large  arrays,  but  these 
are  not  of  interest  to  us. 
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An  undesirable  effect  of  unequal  spacing  is  an  Increase  in  sidelobe 
level.  Figure  3.1  shows  the  spacing  of  the  four-element,  12. IX  QUICK  LOOK 
array.  Figure  3.2  shows  the  corresponding  array  factor  (uniform  weighting). 
While  there  are  no  perfect  grating  lobes,  there  are  several  extremely  high 
sidelobes.  The  performance  bounds  obtained  in  Section  II  show  that  the 
sidelobes  of  the  array  factor  are  very  important  in  determining  the  direction- 
finding  capability  of  an  array. 

C.  Linear  Dependence  of  the  Direction  Vectors  and  Array  Ambiguities 

The  study  of  the  Ziv-Zakai  bound  leads  one  to  the  conclusion  that  when 
direction  vectors  are  nearly  linearly  dependent,  estimation  accuracy 
suffers.  In  this  section,  the  consequences  of  linear  dependence  are 
discussed,  and  a  way  of  testing  an  array  for  this  condition  is  proposed. 

Suppose  there  exists  some  set  of  K  direction  vectors  which  are  linearly 
dependent,  i.e.,  there  exist  complex  constants  such  that 

K 

l  ouyCu.)  ■  0 

k-i  K  K 

This  means  that  the  array  is  completely  blind  to  this  particular  distribution 
of  sources.  Of  course,  the  proper  combination  of  complex  amplitudes  is 
required  in  order  to  get  complete  cancellation;  if  the  sources  are  Incoherent, 
this  situation  would  not  persist  over  multiple  snapshots.  However,  since  any 
one  of  the  signals  can  be  expressed  as  a  linear  combination  of  K-l  others,  a 
distribution  of  K  sources  with  arbitrary  amplitudes  cannot  be  distinguished 
from  any  one  of  K  distributions  of  K-l  sources. 

A  simple  example  of  this  is  a  uniform  linear  array  with  element  spacing 
X .  Such  an  array  has  grating  lobes  at  u"±l .  A  signal  arriving  from  direction 
u  is  indistinguishable  from  one  arriving  from  directions  u±l  (one  of  which  is 
in  "visible  space").  These  signals  are  linearly  dependent.  An  estimator 
faced  with  this  situation  (and  no  a  priori  information)  can  do  no  better  than 
to  choose  one  of  the  two  possible  directions  at  random.  The  result  is  a  large 
mean  squared  estimation  error. 
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Fig.  3.1.  Quick-look  array 
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It  can  be  shown  that  for  equally  spaced  linear  arrays  without  grating 
lobes,  every  set  of  distinct  direction  vectors  Is  linearly  Independent.  This 
result  does  not  extend  to  nonequally  spaced  arrays.  The  following  argument 
shows  why. 

Consider  an  M-element  equally  spaced  linear  array  illuminated  by  L<M 

sources  equally  spaced  in  angle  in  the  unit  circle,  l.e., 
X  1 

u!  "  I  l'  *  "  1,  •  T*ie  re8Pon8e  ®rf*y  element  m  to  source  £  is 

j2ir^r 

x  (1)  -  o.  e  L 
in  jc 


where  a£  Is  the  complex  amplitude  of  source  £ , 
all  sources  is 


The  response  of  element  m  to 


L  L  .-i rt 

l  m  l  «»eJ  L 


£-1 


£-1 


2ir 


Assume  *  1  .  Since  e^  L  is  an  Lth  root  of  unity,  it  follows  that 

,  0  m*0  mod  L 

Xm  L  m-0  mod  L 


If  the  elements  whose  Indices  are  0  mod  L  are  removed  from  the  array,  the 
remaining  elements  all  have  zero  response,  and  the  source  distribution  cannot 
be  detected.  This  Is  just  another  way  of  saying  that  there  exists  a  linear 
combination  of  direction  vectors  (sources)  which  is  identically  zero.  The 
elements  to  be  removed  can  be  varied  by  adjusting  the  relative  phases  of  the 
sources. 

In  order  to  cause  problems,  the  linearly  dependent  source  directions  must 
all  be  in  "visible  space".  The  L  sources  span  an  interval  in  u  of  ^  ; 

this  length  must  be  less  than  2  if  all  sources  are  to  be  in  visible  space. 
Thus  If  the  element  spacing  of  the  array  satisfies  d  >  y  (1-|-)  ambiguity 
problems  will  occur. 
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As  a  staple  example ,  consider  a  five  element  array  with  spacing  d“4-  , 

i  * 

and  three  sources  at  u-0,  ±  Wj  .  The  phase  shift  per  eleaent  serosa  the  array 

is  0°,  ±120°  respectively  for  the  three  sources.  Figure  3.3  shows  how  the 

complex  amplitudes  of  the  sources  add  up  at  each  array  eleaent  when  they  are 

In  phase  at  the  center  eleaent.  When  the  center  eleaent  Is  missing,  the  array 

cannot  distinguish  between  one  signal  at  broadside  and  two  signals  of  the  saae 
2 

amplitude  at  u~±  both  having  180s  relative  phase. 

More  generally,  since  for  Uj"0,  U2*  e  ^1r^,  u3“  ^ 

+  Z2  +  V3  -  0 

any  combination  of  three  sources  from  these  three  directions 
avj  +  bv^  +  c v^ 

cannot  be  distinguished  from  pairs  of  sources 

(a-c)vj  +  (b-c)v2 
or  (a-b)jm  +•  (c-bjv^ 

or  (b-a)^  +  (c-a)v^. 

If  no  information  about  the  slgnal-in-space  covariance  matrix  P  Is 
available,  these  ambiguities  are  fundamental;  there  Is  no  way  to  determine  the 
true  source  directions.  However,  If  the  signals  are  known  a  priori  to  be 
Incoherent  (diagonal  P  aatrix)  the  true  source  directions  can  be  determined. 
The  reason  Is  that  only  one  of  the  possible  configurations  of  sources  will 
yield  e  diagonal  P  aatrix. 

An  algorithm  Incorporating  this  feature  is  currently  under  Investigation. 


Testing  for  Linear  Dependence 

A  set  of  vectors  _v^ ,  Vy ,  .  .vy  is  linearly  dependent  if  and  only  if  the 
matrix  G-V^V  I Vj  Vg  ...  hjj  Vj  ...  VjJ  is  singular,  l.e.,  if  | V**V |-»0. 
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Fig.  3.3.  Example  of  linearly  dependent  direction  vectors:  expurgated 
linear  array. 


I 

» 


i 

i 


The  test  for  linear  dependence  is  simply  to  evaluate,  the  determinant  as  a 
function  of  the  direction  parameters  u^ and  look  for  zeros.  Since  0< | V^V | < 1  *  * 

an  equivalent  procedure  is  to  look  for  peaks  of  1-  j  V^V  | »  For  arrays  of 
Isotropic  elements  and  k-2  signals,  this  function  is  simply  the  magnitude  of 
the  normalized  array  pattern,  and  the  peaks  represent  grating  lobes. 

For  k-3  signals  on  an  isotropic  array,  we  have 


1  fUj-Uj)  f(u3— Uj) 

f*(u2"ul)  ^  1  f(u3-u2) 

f  (Uj-Uj)  f  (u3“U2)  1 


where  f(u)  ■  rr  )  e^  1TXmu  is  the  array  pattern.  This  is  a  function  of 
m*l 

two  independent  variables,  x~U2-ui»  and  Thus 

G(x,y)  -  |f (x) | 2+| f (y) |2+| f fr+y) |2-2Re{  f (x)f (y)f*(x+y)} 

Without  loss  of  generality,  we  may  assume  Uj  <u2<u3,  which  implies  x,  y  > 
0.  We  are  only  interested  in  source  distributions  which  are  entirely  in 
"visible  space",  i.e.,  u3  -  Uj  ■  x+y<2.  Thus  the  region  of  Interest  becomes 
0<y<  2-x),  CK x<2 .  Furthermore,  it  follows  immediately  from  the  relation  f(-u) 
■  f*(u)  that  G(y,x)  •  G(x,y).  We  plot  it  over  the  entire  region  for  esthetic 
reasons. 

Figure  3.4  shows  a  contour  plot  of  the  function  G(x,y)  over  this  region 
for  the  expurgated  linear  array  of  Fig.  3.3  with  element  spacing  1/2.  The 
contours  go  from  0.5  to  1  in  steps  of  0.05.  For  any  linear  array,  the 
direction  vectors  are  dependent  when  two  of  the  directions  are  equal.  This 
creates  ridges  of  ambiguity  along  the  lines  u£~uia0  and  the  left  side 


This  assumes  an  ideal  array  whose  direction  vectors  satisfy  v^(u)v(u)»l. 
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and  I  -.tom  of  the  triangle.  Also,  for  linear  arrays  (expurgated  or  not)  whose 
element  spaclngs  are  multiples  of  A/2,  there  is  a  ridge  of  ambiguity  along  the 
hypotenuse  of  the  triangle,  which  corresponds  to  The  plot  shows 

there  is  also  ambiguity  at  x-y-2/3  which  corresponds  to  the  source 
distribution  previously  discussed. 

The  corresponding  plot  for  the  unexpurgated  array  is  shown  in  Fig.  3.5. 
As  expected,  the  ambiguity  at  x“y*2/3  has  disappeared. 

Somewhat  more  interesting  results  are  obtained  for  the  QUICK  LOOK 
array.  Figure  3.6  is  an  ambiguity  plot  for  this  array.  The  contour  levels 
begin  at  0.7  and  go  upward  in  steps  of  0.05.  The  contours  of  high  ambiguity 
are  predominantly  straight  lines,  either  horizontal  (constant  U3  -  U2>, 
vertical  (constant  U2  -  Uj),  or  sloping  downward  to  the  right  (constant  U3  - 
uj ) .  Close  examination  shows  that  the  peaks  (excluding  those  along  the  line 
u3  -  u2  ■  0)  never  reach  unity.  The  largest  peak,  located  at  ^“Uj  *  u3“u2  “ 
0.33038,  has  value  0.99990. 

The  nature  of  Fig.  3.6  can  be  largely  explained  on  the  basis  of  the  array 
pattern.  This  pattern  is  shown  in  Fig.  3.7  with  the  locations  and  values  of 
the  9  largest  sidelobes  shown.  Lines  representing  the  sldelobe  separations 
have  been  drawn  on  the  figure.  The  ridges  of  high  ambiguity  follow  these 
lines  and  their  intersections  produce  the  largest  peaks.  For  example,  there 
is  a  peak  value  of  0.997  at  U2~Uj“1.65,  U3”U2*0.33.  This  corresponds  to  Uj»- 
1.65,  U2**0,  U3-O.33;  signals  in  the  main  lobe  and  two  large  sidelobes. 

The  ridges  represent  near  linear  dependence  between  two  direction 
vectors.  At  their  intersections  all  three  vectors  are  nearly  equal.  In 
addition,  there  are  additional  peaks,  such  as  the  one  at  u2~ul“u3-u2"®*^* 
where  a  linear  combination  of  the  three  vectors  in  nearly  zero,  but  no  linear 
combination  of  two  is. 

To  predict  the  location  of  these  near  ambiguities  note  that  the  QUICK 

LOOK  array  spacing  is  nearly  3,  2.5,  6.5A.  This  is  an  expurgated  linear  array 

with  y"  “  0*5  and  element  spacing  6,  5,  13  d.  Suitably  phased  equal 
^  d  1 

amplitude  sources  at  —  u-0 ,  ±  produce  outputs  on  the  array  at  element 
numbers  3k+l,  none  of  which  are  present;  consequently,  these  direction  vectors 
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00*2 


Fig.  3.5.  G(x,y)  for  uniform  linear  array. 


are  linearly  dependent.  Although  certainly  not  obvious  at  the  outset,  an 
ambiguity  also  exists  for  j-  u«0,  ±  viz. 


The  array  is  also  nearly  3,  2.4,  6.6X,  which  is  an  expurgated  linear 

array  with  y  **0.6  and  element  spacing  5,  4,  lid.  Sources  at  u*=0, 

XI  \  *  2 

•2-  (-=-+£)  —  (—  +m)  for  any  integers  i,  m  will  produce  zero  input.  An 

d  3  d  3 

example  that  falls  within  the  region  plotted  is  ^-u^uyn^  *  5/9. 

A  third  approximation  to  the  array  is  3.3,  2.2,  6.6X,  with  d/X  *  1.1  and 
spacing  3,  2,  6d.  This  produces  ambiguities  at  U2_uj“U3-U2*0.303,  0.606, 
0.909  at  u2“Uj*1.212,  u2-Ui*1 .212,  u3-u2“0,303,  and  perhaps  elsewhere. 

Since  these  are  only  approximations  to  the  QUICK  LOOK  array,  its 
ambiguities  are  not  perfect.  However,  all  produce  ambiguity  peaks  which 
exceed  0.95. 

The  fact  that  the  ambiguities  of  the  array  are  not  perfect  means  that  at 
sufficiently  high  S/N  ratio,  the  correct  source  distribution  can  be 
determined.  However,  in  near  ambiguous  situations,  the  required  S/N  can  be 
very  large.  This  can  be  shown  via  the  Ztv-Zakai  bound  for  two  signals  (see 
Appendix  VII). 

Array  ambiguities  involving  more  than  three  directions  also  exist.  For  N 
directions,  the  quantity  to  be  examined  is  a  function  of  N-l  variables,  and 
finding  the  peaks  is  a  much  more  laborious  process.  For  this  reason,  we  con¬ 
fine  our  attention  to  ambiguities  involving  only  three  directions. 

This  study  of  linear  dependence  of  direction  vectors  has  led  to  the 
following  conclusions  for  unequally  spaced  arrays. 

1.  Ambiguous  source  distributions  involving  combinations  of 
two  or  more  sources  exist  for  expurgated  linear  arrays. 

2.  For  more  general  unequally  spaced  arrays,  these 
ambiguities  are  usually  not  perfect.  However,  the 
threshold  Bignal-to-nolse  ratio  necessary  to  resolve  them 
can  be  quite  large. 


IV.  QUICK  LOOK  Performance  Simulation 


A.  One  Emitter 

QUICK  LOOK  performance  was  evaluated  first  for  the  case  of  a  single 
emitter.  Theory  predicts  and  simulations  confirm  that  at  high  S/N  ratio,  the 
maximum  likelihood  estimate  of  u  ■  sind  Is  asymptotically  efficient,  l.e., 
approaches  the  Cramer-Rao  bound,  and  thus  is  the  best  possible  (unbiased) 
estimate.  The  ML  estimate  is  just  the  value  of  u  for  which  the  steered  array 
output 


x 

M  -j  2it“  u 

|  l  z  d  X  I 
1  L,  m  ' 

m**l 


is  a  maximum.  This  is  referred  to  here  as  beamsum  processing. 

Current  QUICK  LOOK  processing,  which  assumes  only  one  emitter  is  present, 
uses  only  phase  differences  between  three  pairs  of  elements.  This  processing 
was  also  simulated.  A  correction  to  the  algorithm,  detailed  in  Appendix  D, 
was  found  to  be  necessary  in  order  to  get  correct  results  for  emitters  near 
end  fire. 

Figure  4.1  shows  simulation  results  for  beamsum  processing  for  emitters 
at  u  -  0,  0.5,  .855,  .995  (9  *  0,  30°,  60°,  84°).  The  RMS  error  in 

beamwidths*  is  plotted  vs.  array  signal-to-nolse  ratio  for  each  case.  Also 
shown  for  comparison  is  the  single  emitter  Cramer-Rao  bound.  The  threshold 
SNR,  where  the  error  Increases  abruptly,  depends  on  the  emitter  direction. 
This  Is  due  to  the  nature  of  the  array  pattern  of  the  QUICK  LOOK  array  (Fig. 
3.7),  as  we  now  show. 


*  Hzi  A 

Beamwldth  is  defined  to  be  ^  L  where  M  is  the  number  of  array  elements  and 
L  is  the  array  length. 
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Beamsum  processing  scans  the  beam  over  all  of  "visible  space,"  i.e.,  from 
u  ■  -1  to  u  *  1,  and  takes  the  direction  of  maximum  response  as  the  emitter 
direction.  When  the  emitter  is  at  broadside  (0°),  the  result  of  the  scan  is  a 
noisy  version  of  the  array  pattern  from  u  ■  -1  to  u  ■  1  (the  array  pattern  is 
symmetric  abut  u  *  0).  The  largest  sidelobe  in  this  Interval  is  0.88535  at  u 
■  10.33.  When  the  SNR  drops  sufficiently  so  that  these  three  peaks  can  be 
confused,  the  direction-finding  error  starts  to  increase  abruptly.  When  the 
emitter  is  at  30°,  the  region  of  Interest  is  u  ■  -1.5  to  u  *  0.5.  Since  the 
largest  sidelobe  is  unchanged,  the  threshold  SNR  remains  nearly  the  same. 
With  an  emitter  at  60°,  however,  the  region  of  interest  is  u  ■  -1.866  to  u  * 
0.134,  and  the  largest  sidelobe  is  0.95467  at  u  *  -0.784.  The  threshold  SNR 
should,  therefore.  Increase  by  about  (1-0.88636)^/(1-0.95467)^  or  about  8 
dB.  An  emitter  at  84°  brings  an  even  larger  sidelobe  with  height  0.96962  into 
the  visible  region.  The  resulting  change  in  threshold  SNR  should  be  about  3.5 
dB.  These  predictions  are  in  reasonable  agreement  with  the  simulation 
results. 

The  ZZ  worst-case  bound  (2.33)  does  a  reasonable  job  of  predicting  the 
threshold  SNR.  Figure  4.2  shows  again  the  simulation  results  for  beamsum 
processing  with  a  single  emitter  at  0°  or  at  60®,  together  with  the 
appropriate  ZZ  bound.  For  example,  when  the  emitter  is  at  60°,  the 
appropriate  ZZ  bound  is  the  worst-case  value  taken  over  u  -  -1.866  to  u  ■ 
0.134. 

The  QUICK  LOOK  processing  algorithm  is  not  a  pure  estimator  but  combines 
estimation  and  detection  (see  Appendix  VIII).  When  phase  ambiguity  resolution 
falls,  the  algorithm  rejects  the  data  and  makes  no  estimation  of  direction. 
At  low  S/N  (<  10  dB),  only  about  20%  of  the  trials  result  in  estimates.  Such 
"edited"  estimates  are  not  subject  to  the  performance  bounds  discussed 
previously,  except  in  the  limit  of  large  S/N. 

For  this  reason,  the  QUICK  LOOK  algorithm  was  modified  for  purposes  of 
comparison  with  other  pure  estimators  and  bounds.  When  ambiguity  resolution 
fails,  the  modified  algorithm  chooses  a  random  number  between  -1  and  1  for  the 
emitter  direction  u.  Simulation  results  for  this  modified  algorithm  are  shown 
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Fig.  4.2.  Comparison  of  beam-sum  processing  with  Ziv-Zakai  bounds 


,g.  4.3.  Single  emitter  simulation  results  for  QUICK-LOOK  array  with 
JICK-LOOK  (phase  only)  processing. 


in  Fig.  4.3  together  with  the  CR  bound.  At  high  SNR,  the  algorithm  achieves 
the  CR  bound;  its  threshold  SNR  is  somewhat  higher,  particularly  for  emitters 
near  broadside. 

The  most  important  conclusions  to  be  drawn  from  the  simulations  are: 

(1)  Beamsum  processing  offers  little  improvement  in  the 
single-emitter  case.  The  (modified)  QUICK  LOOK  algorithm 
is  essentially  optimum. 

(2)  Ziv-Zakai  bounds  can  be  used  to  determine  a  bound  on 
threshold  S/N  for  estimating  the  direction  of  arrival  of 
a  single  signal. 

B.  Two  Emitters 

4k 

Two  emitter  CR  bounds  on  the  standard  deviation  of  the  angular  error  are 
shown  as  functions  of  emitter  separation  in  Fig.  4.4.  Plotted  are  the  bounds 
for  a  Gaussian  signal  model  and  for  a  deterministic  signal  model  with  best  and 
worst  case  phase  relationships.  The  received  signals  have  equal  power  (20  dB 
array  S/N)  and  100  observations  are  available.  In  the  deterministic  model, 
the  phase  difference  between  the  signals  is  a  critical  parameter.  The 
difference  in  performance  between  best  and  worst  case  phase  conditions  is 
often  appreciable.  However,  it  was  pointed  out  in  section  II.l.E  that  when 
multiple  observations  are  available  and  the  signal  phases  are  random,  the 
correlation  between  the  directional  errors  will  be  small  and  the  deterministic 
and  Gaussian  signal  models  yield  similar  bounds  on  performance.  For  these 
reasons,  the  Gaussian  signal  model  was  used  in  the  simulation. 

In  the  single  emitter  simulation,  the  maximum  likelihood  (ML)  algorithm 
(beamsum)  was  used.  For  the  two-emitter  simulation,  the  MUSIC  algorithm  [10] 
was  used. 

Figure  4.3  shows  the  results  for  signal  separations  greater  than  1 
beamwidth.  Each  point  represents  the  standard  deviation  in  the  angle  error 
for  a  particular  emitter  after  100  Monte  Carlo  trials.  Two  emitters  of  equal 


*More  precisely,  the  standard  deviation  of  the  error  in  u  ■  sin  B. 
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.5.  Two-emitter  simulation  results  and  CR  bounds  for  QUICK-LOOK 
(emitter  separation  >  beamwidth) • 


power  comprised  the  simulated  signal,  so  two  points  were  produced  for  each  set 
of  100  trials.  Signal  separations  of  0,1  to  0.5  In  u  -  sln6,  and  S/N  ratios 
of  10  to  40  dB  were  simulated.  Shorn  also  are  the  CR  bounds  (Gaussian  signal 
model).  The  agreement  is  very  good,  indicating  that  the  MUSIC  algorithm 
achieves  nearly  optimal  performance  for  signal  separations  in  excess  of  1 
beam width. 

Fig.  4.6  shows  results  for  signal  separations  of  0.05,  0.1,  0.2,  0.4,  and 
0.8  beamwidths.  A  breakdown  in  estimation  accuracy  occurs  as  the  emitters  come 
closer  together  in  angle.  The  mechanism  is  that  the  emitters  are  no  longer 
resolved;  they  produce  only  one  peak  in  the  spectrum.  Consequently  a  spurious 
peak  is  accepted  and  its  angle  estimate  assigned  to  one  of  the  two  emitters, 
causing  a  large  RMS  error.  A  similar  effect  was  noted  in  (3]  where  a  two 
emitter  simulation  was  run  using  the  MUSIC  algorithm  and  a  uniform  linear 
array. 

In  the  absence  of  a  Ziv-Zakai  bound  for  the  two  emitter  case,  it  is  not 
possible  to  say  on  the  basis  of  the  present  results  whether  or  not  the 
resolution  performance  of  QUICK  LOOK  can  be  Improved  by  the  choice  of  a 
different  high  resolution  algorithm.  Simulations  using  other  algorithms  are 
now  being  conducted  to  resolve  this  question. 


V .  Summary 

The  problem  of  evaluating  the  emitter  location  capability  of  arrayB  has 
been  addressed  by  studying  Cramer-Rao  and  Zlv-Zakal  bounds  on  angle  estimation 
accuracy.  These  bounds  are  useful  for  two  reasons: 

1.  They  provide  a  standard  against  which  angle  estimation  algorithms  can 
be  judged. 

2.  They  permit  comparisons  of  various  array  geometries  which  are 
algor i thm-i ndepe nde  nt . 

Although  not  leading  directly  to  an  analytical  theory  of  optimum  array 
geometries,  the  study  clarified  the  notion  of  array  ambiguities  in 
nonunlformly  spaced  linear  arrays  and  revealed  a  method  for  finding  them.  A 
power  estimation  technique  for  resolving  these  ambiguities  when  the  signals 
are  Incoherent  is  now  being  tested. 

The  study  also  showed  the  connection  between  high  array  sldelobes  and  the 
threshold  signal-to-nolae  ratio  required  to  locate  one  emitter.  When  one 
sidelobe  Is  dominant,  threshold  S/N  is  directly  related  to  the  height  of  the 
sidelobe;  this  becomes  less  true  if  multiple  sldelobes  of  roughly  equal  height 
are  present.  Nevertheless,  a  useful  rule  of  thumb  for  reducing  the  required 
S/N  is  to  minimize  the  maximum  sidelobe  of  the  array. 

Special  attention  was  devoted  to  simulating  the  performance  of  the  QUICK 
LOOK  array.  For  a  single  emitter,  QUICK  LOOK  processing  was  found  to  be 
essentially  optimum.  For  two  emitters,  the  lack  of  an  appropriate  Zlv-Zakal 
bound  makes  it  difficult  to  say  whether  the  breakdown  in  emitter  resolution 
capability  is  a  property  of  the  algorithm  employed  or  a  fundamental  limitation 
caused  by  the  array  geometry.  This  question  is  presently  being  resolved  via 
analysis  and  simulation. 


APPENDIX  I 
IMPROVED  ZZ  BOUND 


Wax  and  Zlv  [9]  have  Improved  the  basic  Zlv-Zakai  bound  by  a  factor  of 
two.  Their  derivation  assumes  that  an  estlaate  x  of  soae  parameter  xe[-A,  A] 

A  A  A 

takes  values  in  the  finite,  symmetric  Interval  l -A,  A]  where  A  <  A.  The 
result  is  then 


j  le2(a)  +  e2(-a)j  >  2  a2  Pg,(-a,  «) 

a  bound  on  the  average  mean  squared  error  when  the  true  values  are  x-a  and  x— 
a.  This  is  not  as  specialized  a  result  as  it  appears  to  be  at  first  glance. 

Given  any  parameter  x  on  a  finite  Interval  [A,  B]  and  any  two  values 

X1  +  x2 

Xj,  %2  e[A,  B] ,  the  linear  transformation  y  »  x - ^ - -  x  -  x  maps  xj 

x2  ”  X1 

into  -a,  *2  into  a»  where  a  * - j - .  The  new  parameter  y  takes  values  on 

the  non-symmetric  Interval  (A  -  xQ,  B-xQ].  We  also  use  the  transformation  to 
generate  a  new  estimate  y  taking  values  in  this  same  Interval.  The  mean 
squared  error  of  the  estimate  is  preserved  under  the  transformation. 

A 

Now  define  a  truncated  estimate  y£  on  the  largest  symmetric  Interval 
contained  in  the  mapping  of  [A,  B],  namely  [-Y,  Y],  where 


Y  -  min{|A  -  xj,  jB-xJ) 


A  A 

If  y  lies  outside  of  this  Interval,  y  ■  AY  as  appropriate.  The  mean  squared 

»  A 

error  of  this  truncated  estimate  will  be  smaller  than  that  of  y  when  the  true 
values  are  Xj,  x2.  Thus, 


Y  l*2-  (x^)  +  e2,  (x2)j  »  i  [e2, (a)  +  e2  (-a)J 
*  x  y  y 


(x  • X  )  ^ 

>  4  [e2A  (a)  +  e2A  <— a) |  >  2a2P_  (-a,  a)  ■ — ^—5^ —  Ps  (x,,x,) 


2  ‘EA'“1*  2' 

x 


This  shows  that  the  result  holds  for  any  finite  Interval,  symmetric  or  not 
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APPENDIX  II 

PROBABILITY  OP  ERROR  IN  CHOOSING  BETWEEN  SIGNALS  OF  KNOWN  COMPLEX  AMPLITUDE 

The  contending  hypotheses  are: 

Hji  £  ■  P  sCuj)  +  £ 

H2:  —  m  p  +  — 

where  p  *  Ae^  is  the  known  complex  amplitude, 

I  s^<u± )  | 2  -  1  , 

and  £  is  a  complex  Gaussian  noise  vector  with  zero  mean  and 
E  £  £T«  0,  E  £  £H-  A. 

A  decision  rule  is  a  division  of  the  complex  M-dlmenslonal  space  a  of  all 
possible  observation  vectors  z_  into  two  parts,  aj  and  a2.  If  the  observation 
falls  in  ap  choose  Hjj  otherwise,  H2«  The  a  priori  probabilities  of  the  two 
hypotheses  are  Pj  and  P2- 

The  probability  of  error  for  this  decision  rule  is 

PE  -  P2Pr{z  e  a,|H2}  +  PjPrfz  e  a^Hj 

m  p2  /  p(£|H2)d£  +  Pj  /  p(a|H1)d£ 

al  a’al 

-  p  1  +  J  [P2p(£|H2)  -  PjpUlHjJjdz 
“l 

To  minimize  Pg,  assign  to  a^  if  and  only  if 
PjpCslHj)  >  P2p(z|H2) 
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An  equivalent  and  more  convenient  fora  la 


In  pUjHj)  -  In  p(rjHj)  >  In  ~- 

This  Is  the  well-known  likelihood  ratio  test.  For  the  problea  at  hand, 

"  |  *  -  pjstu^l2 

p(*|H.)  -  — -  e  1-1,2 

-  1  *M|A| 


and  the  likelihood  ratio  test  becomes 


x  »  2Re  |  p£  A  (£i“  s^}  +  A  (s£A  -  a^A  Sj 


>l 

y  ta~-  \ 

k  Pl 


x  la  a  real  Gaussian  randoa  variable  with  mean 
ax  -  E(x|Hj)  -  A* (•*-•* l^-1 

under  Hj  and 

-  E(x|H2)  -  -Bj 

under  H2«  The  variance  of  x,  which  is  the  saae  under  either  hypothesis,  is 
given  by 

a2  -  E(4Re2{peHA  l(£.  -  s.)}) 


The  probability  of  error  is 


P  «P2Pr{x>  X|H2}  +  P  jPrfx  <  X  |  HA } 
G 


2a 


_  (x  +  Bj) 


-  -V (x  ■  mir 


X  2a2 


- -  lP  2  / 


dx  +  ?  i  J  e 


dxj 


/  2ircf 


X-m, 


_L_  tp2  /  e  dx  +  Pj  / 


x 

2 


dxj 


/  2ir 


X+m, 


m,  +  X  “i  ”  < 

p2ecfc*("^ - )  +  P1  > 


m. 


To  get  the  tightest  booed,  «e  oexi.lte  this  expresses  vrtth  te.pett  to  l 


Since 


P2  1  '  P1  X 

L  -  - ■  e 


pi  pi 


we  have 


P1  -  1/(1  +  e  ) 


-X  ■ 


P2  -  1/(1  +  e  ) 


m,  +  X  ,  .  ,  “1 

'  7T7>  -V-’ +  TTX  "£c*(— > 


PE  (X) 
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The  derivative  is 


dP  .  ■  +  X  n  -  X 

dT"  “  TrTTToshxT  ^erfc*("“5  ')  ~  erfcJ~  )] 


+  7=— l— T 

/  2 wo  1  +  eA 


(■j-  x)‘ 

2 a2 


(«!+  xr 

2a  2 


1  +  e 


-X 


When  the  expressions  for  nij  and  ox  are  substituted  into  this  expression, the 
second  tern  vanishes  and  the  derivative  becomes 


dPg  i  m,  +  X  m,-  X 

dX  “  2(l +  coshX 


j  lerfcJ“ - )  "  erfc*(~ — )) 


This  is  obviously  zero  when  X  -  0;  furthermore,  it  is  clearly  always  negative 
for  X  >  0  and  positive  for  X  <  0.  Thus,  Pg(X)  has  a  single  maximum  at  X  ■  0, 
(equal  a  priori  probabilities)  and 


yx  PE(X)  -  erfc*  ^  -  erfc*  (--) 


(II. 1) 


where  d  ■  |sj  -  sj  f  A_1  is  the  distance  between  the  signals. 

For  aultlple  snapshots  n-1,  N  the  likelihood  ratio  test  becomes 

i2/  H.-*_  H.-l 


*  "n^l  l2RetpiAi  A  ^21“  22  ^  +  |Pnl  (*2^  22  “  2lA  21^1 


If  the  total  signal  power  is  defined  to  be 


then  the  previous  equations  for  the  means  and  variance  of  x  still  apply  and 
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the  result  (II. I)  Is  still  valid.  The  only  effect  of  Multiple  snapshots  on 
the  bound  is  to  increase  the  signal-to-noise  ratio. 


APPENDIX  III 


ERROR  PROBABILITY  FOR  CHOOSING  BETWEEN  SIGNALS  WITH  UNKNOWN  PHASE 


The  contending  hypotheses  ere 

Hj!  _z  -  AeJ*  stuj)  +  £ 

H2:  —  "  Ae^  i.^u2^  +  1 


where  the  ststlstlcs  of  are  as  defined  in  Appendix  II. 

J*i 

We  use  the  result  (II.l)  of  Appendix  II ,  replacing  by  Sj  ,  and 
obtain  immediately 

P_  -  erfc  (^)  (III.l) 

E  /  2 

where  now 


d2  -  |e^l  £1~  e^2  £2|2_j  “ 


.-1 


There  remains  only  the  task  of  maximizing  Pg  over  the  remaining  free  variable 

n 

6.  To  maximize  Pg,  we  minimize  dS 

d2  "  lijl 2_l  +  |s2|2_1-2Re{8jA'1S2ejB} 

A  A 

Obviously,  this  expression  is  minimized  by  making  the  last  term  as  negative  as 
possible. 


dL»-  iii'V  ii2i2-r 

A  A 
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When  N  snapshots  r n,  n-l,N  are  available  the  result  is  again  (III.l) 
where  now 


<»«2  -A 


This  Is  minimized  term  by  term  over  the  nuisance  parameters  8n  to  get 

(Ad)2  .Jk  A2  |«2|2.r  2\±yl*2\) 

A  A 

which  Is  used  in  (£11.1). 
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APPENDIX  IV 
MULTIPLE  SIGNALS 


The  preceding  analysis  applies  virtually  unchanged  to  the  multiple  signal 
problem.  The  hypotheses  are 

j$l  t 

H  :  z  -  Ae  V-Cu^  +  J&2  Aie  — ^ui^  +  £  *  VOO  J>  +  £ 

i*i  t 

H'  :  z  -  Ae  v(u'  )  +  £2  A'  e  v(u'  )  +  £  -  V(u'  )  £'  +  £ 

where  u  is  a  vector  of  angles  of  arrival  and  £  is  a  vector  of  complex 
amplitudes.  Because  of  the  multiplicity  of  parameters,  we  have  changed 
notation;  the  parameter  values  under  the  two  hypotheses  are  now  denoted  by 
unprimed  and  primed  variables.  Also,  to  avoid  degeneracies  caused  by 
reindexing  of  parameters,  we  specify  that 

u^  <  u^  <  •  •  .  <  Uj  and  u'  <  u'  <  .  .  .  <  u'  . 

Proceeding  as  before,  we  conclude  that 


V£’  “  erfc*( 


Ad(j>,  £') 


) 


where 


d2  «  \  lv  2  -  V'  £'  I2.!  -  |vj  +  ^  bA  -  ejV1  -  ^ 
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I 


This  is  to  be  minimised  with  respect  of  y.  J>..  _b' ,  and  uit  uj,  i-2,  .  .  .  I  * 
Define 

W  -  AJ/2  [V  V"  ] 


*1  -  %  ~ 

I  so  that 

i 

| 

i  n  9 

d  ■  |s_j  -  Wa| 

which  is  minimized  by  choosing  a_  ■  W+8j . 

d2  -  |(I  -  WW+)is j | 2  -  |(I  -  WW+)  AJ/2(vj  -  P  vpf2 
•  jwj  -  e-*Ywj|  2-  (wj2  +  |vj|2  -  2Re{Wj  w^  ejY> 
min  d2  -  (wj2  +  | wj  | 2  -  2{wJ  Vj| 

,wH  ^  | 

-  (Iwj!  -  IwjD2  +  2|wj|  1*1  (1  -  -f^TT^f 

The  vectors  Wj,  *  are  the  components  of  A"^2  ,A"^2  v| ,  respectively,  which 

are  orthogonal  to  the  space  spanned  by  the  columns  of  W.  The  distance  d  is  a 
function  of  the  21  directions  of  arrival  u^,  .  .  .  uj,  u{,  .  .  .  u{. 
Maximization  over  these  variables  must  be  done  numerically  and  requires  a  fine 
sampling  grid.  This  appears  to  be  impractical  for  more  than  one  signal. 


^Without  the  ordering  restriction  on  the  u's,  d  can  be  made  zero  by  setting 
u2“ui»  *>2"*^*  »  ui-uj.  bj«l.  This  amounts  to  a  reindexing  of  parameters. 
Rote  also  that  A[«Aj  in  accordance  with  the  discussion  in  Section  II. 2. B. 
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APPENDIX  V 

PROBABILITY  OF  ERROR  IN  CHOOSING  BETWEEN  TWO  GAUSSIAN  SIGNALS 


The  contending  hypotheses  are 


H,  :  a  -  p  v(u.)  4-  f 
1  — n  n-  1  —i 


n  ■  1,  •  •  •  N 


H,  :  Z  -  P +  Ln 
2  — n  i>—  2  — n 

where  Is  the  usual  complex  Gaussian  noise  vector  with  zero  mean  and 
covariance  matrix  AQ,  and  pn  Is  a  complex  Gaussian  random  variable  with 
statistics 


E(pn)  -  E(p*)  -  0,  E|p*|  -  MP 


The  observations  are  therefore  zero  mean  complex  Gaussian  random 
vectors  whose  covariance  matrix  depends  on  which  hypothesis  Is  true. 


*  \  *  \,  +  "Vl 


The  log  likelihood  ratio  Is 


P. P(^l  t  •  o  _i  Pj 

*  "  P^z^.g/Hp-  >V*air 

We  assume  that  AQ  Is  positive  definite  so  that  A_0^2  exist.  Then 
l*4l  -  |A^(I  +  MP  w^A1^  -  |Ao||l  4  MP  w^Jl 
where  w^  ■  A~  ^2  .  The  second  determinant  on  the  right  is  readily  evaluated 


I 


by  enumerating  its  eigenvectors  and  eigenvalues,  with  the  result 


|Ai|  -  (1  +  MPwJw^  |A 


For  brevity,  we  define  -  MP  w^w^. 


The  Inverse  A-J  is  evaluated  using  Schur's  formula. 

w  wH)  A~  l/2 
-i-i ;  A0 


A”1-  rl/2n  -  MP 
Ai  Ao  (1  I+T 


Using  these  results  reduces  the  log  likelihood  ratio  to 


i  -  N  in  (1  +  q2)  -  N  in  (1  +  +  inP1  -  JtnPj 


MP 


nr  i  |  n  , 

T -  }  W,  A  *  2 

+  q,  L  ’—1  O 


1  +  q 


H  ,2 _ WP_  v  I..H  . ~  V?  _  .2 


■  i  mr  v  i  n  . "  n  i  - 

,1  -  r+T"  *  *-2  Ao  v 


Define  the  complex  Gaussian  random  vectors 


/MP~ 


(1  +  q  j )~  ^  Wj 

a  +  q2rv2w« 


rV2, 


and  the  2x2  matrix 


The  likelihood  ratio  test  can  then  be  written  in  the  simple  form 


H  A  T  q2 

)  x  Jx  +  N  in 
L  — n  — n  1  +  q^ 


+  in  —■ 
2 


0 


84 


The  first  term  is  Just  the  difference  of  two  chi-squared  variates  wi'h  2N 
degrees  of  freedoe.  If  the  two  components  of  Xq  were  uncorrelated  (and 
therefore  statistically  Independent),  then  the  chi-squared  variates  would  also 
be  Independent  and  the  probability  of  error  could  be  calculated  as  shown  in 
Appendix  VI.  Since  this  is  not  the  case,  we  seek  a  linear  transformation  T  of 
which  will  give  Independent  components  while  preserving  the  simplicity  of 
the  decision  region. 

We  first  compute  the  covariance  matrix  of  Xq  under  each  hypothesis.  The 
results  are 


E  (x_xtt/H1 ) 


, 1  Va  ■ 

(  1  +  d,*  q12 


Hermltlan 


q2  + 


'12' 


i  +  q. 


,'ii 

C* 

C12 


12 

:22 


E(x_x"/H2) 


ql  +  *q12' 

1  +  qj 

Hermltlan 


1  +  'ih 


1  +  q. 


y*  q 


12 


i 


n 

\  ; 

gl 


where  q^  “  MPw^Vj. 

Define  the  transformed  vector. 


■^n 


^-Tx 

The  covariance  matrix  of  ^  under  is 


■  [i  :] 


Cj 1+2Re(a*c12 )+| a| 2c 
Hermltlan 


22 


*  *  *  *  * 

Cllb  +C12C  ab  +C22*C 

|  b|  2c  11<f2Re(c1  ,bc*)+|  c|  2c 


'12 


22J 
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while  Che  test  statistic  t  ■  }  becomes 


i 


To  preserve  the  form  of  the  test  statistic  as  a  difference  of  chi-squared 
variates,  we  must  choose  b«ca*.  The  test  statistic  then  becomes 

1  H  |"l  0  1 

"TTtf"  [o  -1,1-*] 

We  may  as  well  choose  c*l  so  that 


t  - 


1  -  a 


n  t 

'■  2nJ^n 


(V.l) 


With  these  choices,  the  covariance ‘matrix  under  Hi  becomes 


E(Wi> 


c11+2Re(a*c,,)+|a|2c 


'12J 


*  ,  2k  11  c22 
22  C12U  * 

C12 


a  + 


12 


'12 


Hermitlan 


1 2  * 

| a |  cu  +  2Re(a  c12)  +  c22 


The  off-diagonal  elements  are  made  zero  by  choosing 


a 


(V.2) 


where 

2  cn  2  l(l+q1)(l+q2)]1/2q12 

Cll+C22  q  j+  q2+  qjq2+  |q12|2 
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and  Che  covariance  matrix  becomes 


E<WV 


1+<)i  u,  *  ?  q.+lq.J 

’i+2<i^)2r'<*’i2>+  1*1 


1*1  V24^7)1/2|i*<»*''12)  +  -;-7rf' 


(V.3) 


A  similar  exercise  under  H2  leads  to  the  same  expression  for  X ,  shoving 
that  the  transformation  works  regardless  of  which  hypothesis  is  true.  The 
covariance  under  H2  is 


E(WV  - 


+  2(l^)1/2Re(a  q12)  + 


qI  + 


'12' 


1  +  q 


1  +  q 


q2  +  2(1  -Fq 


i)V2 


Re(a  q12>  +  I  a 


2  ql  +  ^q12^ 


1  +  q 


1  -» 


(V.4) 


We  have  not  yet  chosen  a  sign  in  the  expression  (V.2)  for  a.  Since  any 
covariance  matrix  is  nonnegative  definite,  we  have  | Cj2 1 2  £.  cnc22*  and  80 


2|c 


-  M  -i-r-T 


12' 


>  1  - 


11  22 


2/c11C22 
Cll+  C22 


(/cll-  /c22)i 
cll+  c22 


>  0 


so  that  |A|  <_  1.  From  a  right  triangle  with  hypotenuse  1  and  side  |\|  we  have 
\X\  >  1  -^1  -  |X? 

which  shows  that  |a|  1  provided  we  choose  Che  negative  sign  in  (V.2).  This 
in  turn  means  that  the  factor  1  -  [a|2  is  nonnegative,  and  so  the  algebraic 
sign  of  the  test  statistic  (V.l)  is  preserved. 


We  can  now  compute  the  probability  of  error  for  the  teat.  It  is 
PE  -  Px  Pr{t  <  6/H^  +  (1-PpPrtt  >  6/H2} 

1+ql  P2 
6  .NAn_  +  ln_ 

Pr(t  <  fi/H^  -  Pr(^  |ynl|2  <  l  |yn2 1 2  +  (1  "  |  a|  2)fi ) 

is  computed  using  the  results  of  Appendix  VI.  The  variances  of  ynl  and  yn2 
are  obtained  from  (V.3)  under  and  from  (V.4)  under  H2. 

The  probability  of  error  is,  of  course,  a  function  of  the  a  priori 
probability  that  is  true.  To  get  the  tightest  Ziv-Zakal  bound  we  must 
choose  P^  to  maximize  Pg.  We  now  show  that  if  the  noise  is  white  (Aq  *  I), 

P1  *  V2 is  the  appropriate  choice. 

2 

When  AQ  ”  a  I,  we  have 

MP  A  ,  .H  ,  . 

qi  "  q2  “  “  "  O'  q12  "  ^(ul>  — (u2} 

cr 

and  the  covariance  matrices  become 


where 


and 


E<Wi>  • » 


l+2Re(a  v1^)  +  |a|2 


—1—2 ' 


1  +  p 


|a|2+2Re(a  _v^2) 


*  H  ,  .  1  +pl^2l 


1  +  p 


(V  .5) 


E<V>2> 


.  .  ,  H  .2 

r.  2  *  1+PlliZ2i 

I  a  |  +  2Re(a  v^Vj)  + - j- +"p -  0 

1  ,  I  H  1 2 

»  *  1  +  p  IZ1Z2 1 

LO  lap'4-  2  Re  (a  v^)  + - l~T~p - 


(V.6) 
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We  observe  that  the  variates  yn^ ,  yn2  simply  change  their  statistical 
identities  under  an  hypothesis  change.  Thus  U  ■  E|yn].|^»  v  “  r  | yn2 ( ^  also 
reverse  roles.  It  follows  that  in  any  conditional  probability  expression 
involving  these  variable,  we  can  change  hypotheses  by  Interchanging  U  and  V, 
e.g. 


Pr(U  >  V  +  6/Hj)  -  Pr(V  >  U  +  6/H2) 


Since  for  white  noise  qj  ■  q2»  it  follows  that 

1-P. 


6  «  In 


1 


It  will  be  convenient  to  minimize  Pg  with  respect  to  6  rather  than  Pj .  Thus 
„  _  1 


1  « 

1+e 


I 


1-P  * 

1  ,  .  S 

1+e 

The  probability  of  error  is 

pE  (6) - y-  Pr{U  <  V  +  bS/Hj}  +  Pr{U  >  V  +  bS/fij} 

l+e°  1+e 

where  for  brevity  we  have  put  b  ■  1  -  |a|^.  We  first  show  that  Pg  is  a 
symmetric  function  of  6. 


?-(-«)  -  — i-r  Pr{U  <  V  -  b6/H.}  +  — Pr{U  >  V  -  bfi/H,} 

&  ,  .  -6  1  i.O  l 

1+e  1+e 

-  — ~  Pr{U  >  V  +  b6/H2)  +  — y  Pr{U  <  V  +  M/Rj}  -  Pg(6  ) 


1+e 


1+e 


This  result,  as  well  as  those  to  follow,  uses  the  fact  that  the 
distribution  function  of  U-V  is  continuous,  so  that  the  distinction  between  <, 
<  can  be  ignored. 


The  derivative  of  Pg  is 


- *-= - -r-  lPr{D  >  V  +  b4/H_}  -  Pr{D  <  V  +  bS/H.}  1 

(1+e®  )(l+e"^)  2  1 


(V.7) 


+  b(— ~  Pr{0 
14« 


V  +  W/H,} - W  Pr{D 

1  1+e-5 


V  +  m/h2>] 


It  la  readily  verified  that  this  expression  vanishes  when  4-0,  and  that 


Pg(0)  -  Pr{0  <  V/Hj) 


Furthermore,  when  5  >  0,  the  first  bracket  In  (V.7)  Is  negative.  He  now  show 
that  the  second  bracket  Is  also  negative. 

The  pdf  of  U  under  Hj  Is 


P(0) 


U 

«*» 

1  0*"1 

TFUT  (2o2)S  ' 


0  >  0 


and  8lallarly  for  V.  Define  for  brevity 


x  _  U  y  .  V  r  _  WV  c  _  b4 

K  K  T’  ^0 

the  variance*  being  those  under  Hj  given  by  (V*5)e 

«  ,  i  1  1 2 

2  2,.,  ,2,  WolS*! 

"o  -  «t  ■  <1_  1*1  >  <l - n; — 


TV 

since  |a|  <  1,  | I  <  1  end  therefore  0  <  r  < 


Note  that 

>  0 

1. 
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Then 


and 


Pr{D 


V  +  M/Hj}  -  /  pv(v)Pu(vfK)dY 


1  e  ~  , 

— y - 5“m  J 

2 oj  <(N-1)!)V  o 


y  /  . 


pr{U 


V  +  M/H2}  -  pr{V  -  U  +  bd/Hj} 

-  /  pJJ(u)pv(ufbfi  )du 


r  >“1*  . 

/  x  (x+c)  a 


l-r 

-  e~C  r  Pr{D  -  V  +  bS/Hj} 

The  second  ten  can  then  be  written  as  constant 
which  la  negative  for  w  >  0,  r<l. 


“D. 


He  have  shown  that  the  derivative  of  Pg  vanishes  at  d-0  and  is  negative 


for  6>0.  It  follows  f row  the  syswetrp  of  Pg  that  its  unique  wsxiwuw  occurs  at 
d-OCPi  -V$. 
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APPENDIX  VI 

PROBABILITY  THAT  THE  DIFFERENCE  OP  TWO  X2  VARIATES  EXCEEDS  A  THRESHOLD 


Lee  a  and  v  be  X2  variates  Kith  2H  degrees  of  freedoe,  having  variances 
(per  degree  of  freedoa)  o2,  o2  respectively.  The  pdf  of  a  Is 


p(u)  - 


I 


N-l  2o 
u  u 

75FI7T* 


u  >  o 


m.i) 


Then  for  8  >  0 


Pr(u  >  v  +  8 )  -  /  p(v)  /  p(u)dudv 
o 


(VI. 2) 


Using  (VI. 1)  and  (VI. 2),  we  obtain 


Pr(u  >  v  +  8  )  - - - - *•  /  vN_1e“v  /  u^e^dudv  (VI.3) 

((N-l)ir  o  rv+c 


°v  8 

where  r  «  — y  ,  c  -  — j  .  Using  ths  result 


2o 

u  u 


*  w_i  N-l  n 

/  u^e^du  -  (N-l)!  I  e 

a  n-0 


(VI. 4) 


In  (VI. 3),  «s  expand  (rv+c)n  In  a  blnoelal  series,  use  (VI. 4)  again,  change  a 
itlon  Index,  and  obtain  eventually  (still  for  8  >  0) 


N-l  .  a  N-l-e  n 

Pr(u  >  v  ♦  8)  -  X  (  .  l)  - (X  |r  •  c)  (VI.5) 

»-0  "  (1+r)  ***  n-0  111 


which  Is  the  desired  result. 

If  8  <0,  (VI.3)  can  he  uaed  to  get  the  desired  result. 


Pr(u  >  v  ♦  8)  -  Pr(v  <  u  -  8)  -  1  -  Pr(v  >  u  -  8) 
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The  last  probability  can  be  cospatad  fros  (¥1.5)  since  -4^0. 
Whan  4  -  0,  this  reduces  to 

N-l  .  a 

Pr(u  >  v>  -  l  (  .  *> - £- 

a-0  "  (1+r) 


which  can  be  put  In  sore  convenient  fora. 

Pr(u  >  v)  - - i_  l  (IHJ“1)  r*  (1+r)1 

(1+r;  1  n-0  “ 


l-l-a 


Expand  the  last  factor  In  a  binoalal  aeries. 

N-l 

(l+r)iN_i  sfO 

Interchange  the  order  of  suaaatlon 


N-l-n 


Pr(u  >  ,) - I  (»*■')  r-  £  m 

1  — »  ■  n-0 


Pr(u  >  r)  - 


Replace  n  by  N-l-n 


R-l 

l  r' 


N-l-n 


N-l-n 


1H«-1  x  ,N-n-l, 


(l+r)2"-1  o-O 


i  rz‘l)  c 

a-0  * 


Nfa-1  j 


.  W-l  n  n 

Pr(u  >  v)  - ktt  I  r  l  ( 

(1+r)2"1  n-0  a-0 

Finally,  using  an  Identity  Involving  binoalal  coefficients  (Feller  [13). 
p.62).  we  obtain 

Pr(u  >  v)  -  Y  i2*1)  (3^)°  (1k>2""l~n 
n— 0 

This  is  Just  the  probability  of  less  than  II  successes  In  2E-1  trials  with 
probability  of  success  (y^r)  .  This  fora  has  been  given  previously  by 
Rally  (14). 


93 


APPENDIX  VII 

LINEAR  DEPENDENCE  AID  THE  ZIV-ZAKAI  BODED 


A  dtnlopant  of  tbo  Zlv-Zakal  bound  on  accuracy  in  estimating  direction 
Q|  in  tbo  prooonco  of  nultlpla  signals  is  given  in  Appendix  IV.  It  is  shown 
that  the  probability  of  error  is 


•rf« *( 


whore,  for  two  signals, 

7  *  2 

dZ  -  |<I  *  Ww  )  8j| 

•  A~  lv(uj)  -  e^vCupl  A  vCuj)  -  e^v  (u^) 
W  -  A*^v(u2)  ,  v(u-  j )  1  ^  lv(u2)»  v(u')] 


The  vector  (I  -  WW+)  ^  is  just  the  conponent  of  jj  orthogonal  to  the  subspace 
spanned  by  the  colunns  of  W.  Suppose  that  v(u^),  v(u2>,  v(u^)  are  linearly 
dependent  direction  vectors*  and  set  uj  ■  u2,  Uj  ■  u^.  Then  ij  Is  a  linear 
coablnatlon  of  the  colunns  of  W  and  so  d«0.  This  leads  to  the  two-point  bound 


19  9  Cu«  “  U«) 

^eZ(Uj,A)  +  eZ(u2,A)j  >  -Z- 

since  -  -y.  He  are  aaaured  that  the  error  at  either  uj  or  U2  (whichever  is 
worse)  exceeds  the  bound,  no  nstter  how  large  the  aignal-to-nolse  ratio  nay 
be. 

Suppose  Instead  that  v(u^),  v(u«),  v(u^)  are  "almost"  linearly  dependant 
in  the  sense  that  the  determinant  ff*f|  is  nearly  sero.  The  component  of  Sj, 


*81nce  v!j  »A  v.  and  A’^Zls  nonsingular,  the  v*  are  linearly  Independent  if 
end  only  if  the  vj  ate. 
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or  equivalently  ,  orthogonal  to  the  apace  spanned  hy  Vj,  has  squared 
»  length 

d  ’  Sf8  1—1  -  *— 2  "  *23 1 

■  Bin  wH  VH  V  w  nbjeet  to  •  1 
where  tP  -  (1,  0,  0).  The  solution  to  this  faalllar  problea  is 

i1 _ ! _  .  j£1l 

.W1.  1^1 

where  ^  “  (v..  Vj)*  The  denoelnator  Is  a  Measure  of  the  linear  dependence  of 
%2  *nd  ^ ,  which  does  not  affect  the  accuracy  of  estluating  uj . 

2 

A  contour  plot  of  1-  d  as  a  function  of  U2  -  uj  and  u^  -  u2  for  the 
QUICK  LOOK  array  is  shown  in  Fig.  VII. 1*.  Since  d2  is  not  ey Metric  in  these 
two  variables,  the  entire  region  0  <  Uj  -  u^  <  2  -(u2  -  u,  )  <  2  eust  be 
plotted.  The  result  is  very  sleilar  to  a  contour  plot  of  |v^|  (Fig.  3.6) 
except  that  the  horizontal  ridges  which  represent  near  linear  dependence  of _v^ 
and  vj,  are  eisslng.  This  is  due  to  the  denoelnator  factor.  If  we  coapute  d2 
for  estlaatlng  U3,  the  vertical  ridges  disappear;  for  u2,  the  slanting 
ridges.  Since  we  are  ultleately  Interested  in  ell  three  accuracies,  the 
Granlan  determinant  |v% |  is  aost  appropriate. 

Since  d  >0,  increasing  S/N  ratio  will  eventually  make  PE  snail. 
Comparing  two  situations,  one  with  minimum  distance  dj,  the  other  with  minimum 
distance  d2,  the  second  will  require  20  log  d2/d}  dB  additional  S/N  ratio  for 
equal  accuracy. 

He  have  shown  that  whan  v2  is  nearly  a  linear  combination  of 
3*2  ,nd  V3  th*  Zlv-Zakal  bound  predicts  large  estimation  errors  in  trying  to 
distinguish  between  a  pair  of  signals  from  directions  uj,  02  and  a  pair  from 
u2,  u3. 


*For  white  noise,  l.e.,  AQ- o2i. 
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Fig.  VII. I.  Contour  plot  of  1  -  d  fflln  for  QUICK-LOOK 


APPENDIX  VIII 

QUICK  LOOK  PROCESSING  MODIFICATION 

Phase  differences  are  Measured  over  baselines  of  2.5,  5,5,  and  12.1 
wavelengths  (at  the  highest  operating  frequency  fh)  .  Without  noise,  an 
eaitter  at  frequency  f  and  direction  u  will  produce  the  results 

♦2.5 

*5.5  1  5.5  tT  u  -  N  (VIII. 1) 

♦l2.1  *  12-‘  fr  “  ‘  1 

where  the  phases  are  expressed  In  cycles  (l.e.,  fractions  of  a  complete 
revolution),  fr  Is  the  frequency  ratio  f/f^  end  L,  M,  N  are  Integers.  From 
these,  we  derive 

♦2.5 +  M  *5.5 +  N  5 

*5.5  +  B  '  *12.I+L  "  11 
or 

11  *2  5  -  5  A5>5  -  5N  -  11  M  -  Ij 

(VIII. 2) 

11  *5.5  -  5  *12.1  *  5L  -  11  «  -  I2 

In  reality,  of  course,  the  quantities  on  the  left-hand  side  will  not  be 
Integers  because  of  noise  and  interference.  However,  If  these  effects  are 
snail,  they  will  be  nearly  Integers,  and  by  rounding  than  to  the  nearest 
integers ,  correct  results  will  be  obtained.  As  the  Interference  be rones 
larger,  the  rounded  values  will  be  Incorrect.  Since  5  and  li  are  relatively 
prlne,  equations  (VIII. 2)  always  have  a  solution,  l.e.,  for  any  Integers  Ij, 
I2,  there  exist  Integers  L,  M,  N  such  that  the  equations  are  aatisfled. 
However,  nost  such  solutions  can  bs  discarded  as  physically  lnposslble.  Since 
-1  <  u  <  1,  the  first  equation  (VIII. 1)  implies 

-2.5  fr  <  ♦  +  M  <  2.5  fr 
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Sloe* 


0  <  ♦  <  1  , 
we  have 

-2.5  f,  -  1  <  -2.5  fr  -  ♦  <  M  <  2.5  fr  -  ♦  <  2.5  ff 
Since  M  le  an  Integer, 

-12.5  £rJ  <  M<  [2.5  fr] 

-[5.5  £fJ  <  N  <  [5.5  fr)  (VIII. 3) 

-[12.1  frl  <  L  <  [12.1  frJ 

where  [x]  denotes  the  largest  Integer  <  x. 

Solutions  which  do  not  satisfy  these  constraints  can  be  discarded. 

The  equations  (VIII. 1)  are  solved  as  follows:  the  first  equation  Is 
solved  for  N  aodulo  11.  Since  11M  ■  0  nod  11,  5Nj  =  Ij  s  -10  Ij  nod  11,  so 

Hj  =  -21 j  nod  11  (VIII. 4) 

A  slallar  procedure  applied  to  the  second  equation  yields 

»2  3  -I2  nod  5  (VIII. 5) 

The  physical  constraints  (VIII. 3)  require  (for  fr  •  1)  that  -6  <  H  <  5. 
Equation  (VIII. 4)  can  be  rewritten  in  the  equivalent  forn 

llj  3  (28  -  21  j)  nod  11  -  6  (¥111.4* > 

where  now  -6  <  H j  <  4.  Slnllarly,  (VIZI. 5)  can  be  written 


98 


n2  s  (6  -  I2)  aod  5  -  1  (?III. 5*) 

where  -l  <  N2  <  3.  It  Is  now  necessary  to  check  whether  or  not  ud  H2 
■etch  within  the  Interval  {-6,  5] .  Figure  VIII.  1  shows  that  watching  (fj  with 
1*2 ,  N2+5,  N2-5  covers  all  poeslbllltlaa  but  one,  naaaly  that  Thus,  If 
no  natch  occurs  after  three  comparisons ,  we  must  cowpare  Nj+11  with  H2+5» 
This  is  the  required  modification  to  the  QUICK  LOOK  processing  algorithm.  A 
block  diagram  of  the  modified  processing  algorithm  le  shown  In  Fig.  VIII. 2. 
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Fig.  VIII.  I.  Matching  Sj  and  »2  In  the  QUICK-LOOK  phase  anblguity 
resolution  algorithm. 


Fig.  VIII. 2.  Modified  QUICK-LOOK  processing  algorithm 
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Cramer-Rao  bound 
Ziv-Zakai  bounds 


emitters 

direction-finding  performance 


signal-to-noise  ratio* 
linear  arrays 
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The  evaluation  of  the  direct  ion -finding  capabilities  of  non-uniform  arrays  Is  approached  via 
algorithm-independent  lower  bound*  on  achievable  angle  estimation  errors.  Two  classes  of 
bounds  are  considered.  The  first,  known  as  the  Cramer-Rao  bound,  applies  only  to  unbiased  esti¬ 
mates.  Compact  analytical  expressions  for  these  bounds  are  developed  which  are  applicable  to 
very  general  direction-finding  problems,  including  an  arbitrary  number  of  emitters. 
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It  is  well-known  that  Cramer-Rao  boundi  are  overly  optimistic  at  low  signal-to-noiae  ratios.  As 
this  ratio  decreases,  a  point  b  reached  at  which  estimation  accuracy  decreases  abruptly.  Another  class 
of  bounds,  known  as  Ziv-Zakai  bounds,  provide  information  about  the  location  of  this  threshold 
point.  Their  study  suggests  that  poor  direction-finding  performance  occurs  in  situations  where  the 
emitter  direction  vectors  are  part  of  a  set  which  is  nearly  linearly  dependent. 

Such  linear  dependence  does  not  occur  in  the  case  of  uniformly  spaced  linear  arrays  (without  get¬ 
ting  lobes).  However,  it  does  occur  when  dements  are  removed  from  such  arrays.  A  systematic  test  b 
developed  to  test  a  given  array  geometry  for  this  condition. 

Finally,  direction-finding  performance  of  the  QUICK  LOOK  array  b  evaluated  via  both  bounds 
and  simulation. 


